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1 Introduction 

1.1 Introduction and Overview 

Orthogonal polynomials on the real line (OPRL) were developed in the nineteenth century and 
orthogonal polynomials on the unit circle (OPUC) were initially developed around 1920 by Szego. 
Their matrix analogues are of much more recent vintage. They were originally developed in the 
MOPUC case indirectly in the study of prediction theory [llSl ^\ M [ml [HSl [TlHl [M] in 
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the period 1940-1960. The connection to OPUC in the scalar case was discovered by Krein [ISlj . 
Much of the theory since is in the electrical engineering literature [36l EZl [Ml [Ml SOI El 1120^ 1121^ 
[na [1231 [203]; see also [8^ [861 [871 [88l [112] . 

The corresponding real line theory (MOPRL) is still more recent: Following early work of Krein 
|133j and Berezan'ski [9j on block Jacobi matrices, mainly as applied to self-adjoint extensions, there 
was a seminal paper of Aptekarev-Nikishin [3] and a flurry of papers since the 1990s [lOl [TTl [T2l 
[ll[ll[ni[I3[2Ql[2ll[2a[29l[Ml[l3l[45l[4l 

[Ml[Ml[65l[Ml[671[Ml[Ml[IIl[Z3l[H[Za[Z6l[27,^ 108^ ilMl 

[TTn l [ml [TT21 [m l [T37 1 [139 1 [Tim [Tim [TM [^ 
[1861 [1981 [2001 [20n [Ml [203!; see also [7]. 

There is very little on the subject in monographs — the more classical ones (e.g., |23 [ l82 [ [93 lll84j ) 
predate most of the subject; see, however, Atkinson [5l Section 6.6]. Ismail |118j has no discussion 
and Simon [1671 I168j has a single section! Because of the use of MOPRL in [M]; we became 
interested in the subject and, in particular, we needed some basic results for that paper which 
we couldn't find in the literature or which, at least, weren't very accessible. Thus, we decided to 
produce this comprehensive review that we hope others will find useful. 

As with the scalar case, the subject breaks into two parts, conveniently called the analytic theory 
(general structure results) and the algebraic theory (the set of non-trivial examples). This survey 
deals entirely with the analytic theory. We note, however, that one of the striking developments in 
recent years has been the discovery that there are rich classes of genuinely new MOPRL, even at 
the classical level of Bochner's theorem; see [Ml ISSl HQl |72l [Ml [Ml Ell EH ES [ESI [161] 
and the forthcoming monograph ^3] for further discussion of this algebraic side. 

In this introduction, we will focus mainly on the MOPRL case. For scalar OPRL, a key issue 
is the passage from measure to monic OPRL, then to normalized OPRL, and finally to Jacobi 
parameters. There are no choices in going from measure to monic OP, P„(x). They are determined 
by 

Pn{x) = + lower order, {x\Pn) =0 j = I, . . . ,n - 1. (1.1) 
However, the basic condition on the orthonormal polynomials, namely, 

{Pn,Pm) = Snm (1-2) 

does not uniquely determine the pn (x) . The standard choice is 

/ \ Pn\X) 
Pn[x) = -TTWir- 
1 1 -^n 1 1 

However, if 9q,0i, . . . are arbitrary real numbers, then 

Pn{x) = (1.3) 

also obey ()1.2p . If the recursion coefficients (aka Jacobi parameters), are defined via 

XPn = an+lPn+l + K+lPn + anPn-1, (1-4) 

then the choice ()1.3p leads to 

bn = bn, 5„ = e^^"a„e-^^"-i. (1.5) 
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The standard choice is, of course, most natural here; for example, if 

Pn{x) = KnX^ + lower order, (1-6) 

then a„ > implies k„ > 0. It would be crazy to make any other choice. 

For MOPRL, these choices are less clear. As we will explain in Section \\.2\ there are now two 
matrix-valued "inner products" formally written as 



{{f.g))R = j f{x)Ufi{x)g{x), (1.7) 
{{f,9))L= [ g{x)dix{x)f{x)\ (1.8) 



where now /i is a matrix- valued measure and ^ denotes the adjoint, and corresponding two sets of 
monic OPRL: -P^(x) and P^{x). The orthonormal polynomials are required to obey 



9^,p^))R = 5nml. (1.9) 

The analogue of (|1.3p is 

p^{x) = P^{x){{P^,P^r'/'ar. (1.10) 

for a unitary (T„. For the immediately following, use to be the choice an = ^- For any such 
choice, we have a recursion relation, 

XPnix) = P^+i(x)Ai+i +p^{x)Bn+l+P^_i{x)An (1.11) 

with the analogue of (|1.5|) (comparing cj,„ = 1 to general o"„) 

Bn = CrliBnCTn A„ = CJ^_;^A„Cr„. (1.12) 

The obvious analogue of the scalar case is to pick an = 'i-, which makes Hn in 

Pni^) = i^nx"' + lower order (1-13) 
obey Kn > 0. Note that (jl.ll|) implies 

Kn = Kn+lAl_^_^ (1.14) 

or, inductively, 

Kn = {Ai...A\)-\ (1.15) 

In general, this choice does not lead to An positive or even Hermitian. Alternatively, one can pick 
an so An is positive. Besides these two "obvious" choices, k„ > or An > 0, there is a third that 
An be lower triangular that, as we will see in Section [1.4^ is natural. Thus, in the study of MOPRL 
one needs to talk about equivalent sets of Pn and of Jacobi parameters, and this is a major theme 
of Chapter [2j Interestingly enough for MOPUC, the commonly picked choice equivalent to An > 
(namely, pn > 0) seems to suffice for applications. So we do not discuss equivalence classes for 
MOPUC. 

Associated to a set of matrix Jacobi parameters is a block Jacobi matrix, that is, a matrix 
which when written in / x / blocks is tridiagonal; see (j2.29p below. 

In Chapter El we discuss the basics of MOPRL while Chapter [3] discusses MOPUC. Chapter [H 
discusses the Szego mapping connection of MOPUC and MOPRL. Finally, Chapter [5] discusses the 
extension of the theory of regular OPs |180j to MOPRL. 

While this is mainly a survey, it does have numerous new results, of which we mention: 
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(a) The clarification of equivalent Jacobi parameters and several new theorems (Theorems 12.81 and 

(b) A new result (Theorem 12.280 on the order of poles or zeros of m(z) in terms of eigenvalues of 
J and the once stripped J^^^ . 

(c) Formulas for the resolvent in the MOPRL (Theorem [229]) and MOPUC (Theorem 13. 24|) cases. 

(d) A theorem on zeros of det($^) (Theorem 13. 7p and eigenvalues of a cutoff CMV matrix (The- 
orem I3.10p . 

(e) A new proof of the Geronimus relations (Theorem 14. 2p . 

(f) Discussion of regular MOPRL (Chapter [5]). 

There are numerous open questions and conjectures in this paper, of which we mention: 

(1) We prove that type 1 and type 3 Jacobi parameters in the Nevai class have An — > 1 but do 
not know if this is true for type 2 and, if so, how to prove it. 

(2) Determine which monic matrix polynomials, can occur as monic MOPUC. We know 
det{^{z)) must have all of its zeros in the unit disk in C, but unlike the scalar case where 
this is sufficient, we do not know necessary and sufficient conditions. 

(3) Generalize Khrushchev theory [l25l [Ml ITOj to MOPUC; see Section [3l3l 

(4) Provide a proof of Geronimus relations for MOPUC that uses the theory of canonical moments 
[33] ; see the discussion at the start of Chapter [H 

(5) Prove Conjecture [5J] extending a result of Stahl-Totik [THO] from OPRL to MOPRL. 

1.2 Matrix- Valued Measures 

Let Ail denote the ring of all / x Z complex-valued matrices; we denote by the Hermitian 
conjugate of a G Aii. (Because of the use of * for Szego dual in the theory of OPUC, we do not 
use it for adjoint.) For a G Mi, we denote by ||a|| its Euclidean norm (i.e., the norm of a as a linear 
operator on with the usual Euclidean norm). Consider the set V of all polynomials in z S C 
with coefficients from Aii. The set V can be considered either as a right or as a left module over 
Aii; clearly, conjugation makes the left and right structures isomorphic. For n = 0,1, . . . , Vn will 
denote those polynomials in V of degree at most n. The set V denotes the set of all polynomials in 
z £ C with coefficients from C'. The standard inner product in is denoted by (•, ■)£i. 

A matrix-valued measure, on M (or C) is the assignment of a positive semi-definite / x I 
matrix ii{X) to every Borel set X which is countably additive. We will usually normalize it by 
requiring 

/i(M) = 1 (1.16) 

(or /u(C) = 1) where 1 is the I x / identity matrix. (We use 1 in general for an identity operator, 
whether in Mi or in the operators on some other Hilbert space, and for the zero operator 
or matrix.) Normally, our measures for MOPRL will have compact support and, of course, our 
measures for MOPUC will be supported on all or part of dD (B is the unit disk in C). 
Associated to any such measures is a scalar measure 

f,tr{X) = TT{f,{X)) (1.17) 

(the trace normalized by Tr(l) = I), /xtr is normalized by /itr(IK) = ^• 

Applying the Radon-Nikodym theorem to the matrix elements of //, we see there is a positive 
semi-definite matrix function Mij{x) so 



dfiij{x) = Mij{x)dfitv{x)- 



(1.18) 
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Clearly, by (fTTTl) . 

Tr(M(x)) = 1 (1.19) 

for d/itr-a.e. x. Conversely, any scalar measure with /Utr(K) = I and positive semi-definite matrix- 
valued function M obeying (jl.lOp define a matrix- valued measure normalized by (jl.l6p . 
Given / x I matrix- valued functions f,g, we define the I x / matrix {{f,g))R by 



{{f,9))R = J f{x)^M{x)g{x)df,Ux), (1.20) 

that is, its (j, k) entry is 

/ fnj{x) Mnm{x)gmk{x)dfitr{x). (1.21) 

nm 

Since f^Mf > 0, we see that 

((/,/))/?. >0. (1.22) 

1/2 

One might be tempted to think of ((/, f))ji as some kind of norm, but that is doubtful. Even if 
is supported at a single point, xq, with M = l"^!, this "norm" is essentially the absolute value of 
A = /(xo), which is known not to obey the triangle inequality! (See \1Q9\ Sect. I.l] for an example.) 
However, if one looks at 

\n = {Tv{{f,f))n)'/', (1.23) 



one does have a norm (or, at least, a semi-norm). Indeed, 

{f,g)R = TT{{f,g))R (1.24) 

is a sesquilinear form which is positive semi-definite, so (jl.23p is the semi-norm corresponding to 
an inner product and, of course, one has a Cauchy-Schwarz inequality 

\Tr{{f,g))R\ < WfUhWR. (1.25) 

We have not specified which /'s and ^'s can be used in (ll.20p . We have in mind mainly 
polynomials in x in the real case and Laurent polynomials in z in the (?D case although, obviously, 
continuous functions are okay. Indeed, it suffices that / (and g) be measurable and obey 



Trifix)f{x))dMx) <^ (1.26) 

for the integrals in (11.2ip to converge. The set of equivalence classes under / ~(7if||/ — 5(||/j = 
defines a Hilbert space, Ti., and {f,g)R is the inner product on this space. 
Instead of (ll.20p . we use the suggestive shorthand 



{{f,g))R = J f{xVdfi{x)g{x). (1.27) 
The use of R here comes from "right" for if a £ Mi, 

{{f,9a))R = {{f,9))Ra, (1.28) 

{{fa,g))R = amf,g))R, (1.29) 
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but, in general, {{f,ag))R is not related to {{f,g))R- 

While (Tr((/, /))_r)^/^ is a natural analogue of the norm in the scalar case, it will sometimes be 
useful to instead consider 

[det((/,/))H]i/2. (1.30) 

Indeed, this is a stronger "norm" in that det > =^ Tr > but not vice- versa. 

When dfi is a "direct sum," that is, each M{x) is diagonal, one can appreciate the difference. 
In that case, dfi = dfii © • • • © d^ii and the MOPRL are direct sums (i.e., diagonal matrices) of 
scalar OPRL 

P^(x, d^x) = Pn{x, dfXi) © • • • © Pn{x, dfll). (1.31) 

Then 

\\Pn\\R = [T.\\Pn{-,dti,)\\l,^,A , (1.32) 
^ 7 = 1 ^ 



while 



(det((P„«,P„«))fi)i/2 = ll\\Pn{;df,j)h2^,^^y (1.33) 



In particular, in terms of extending the theory of regular measures |18Uj . H-P^H]/" is only sensitive to 
max||P„(-, d//j)||]^2^^^^.') while (det((P^, P^f ))/j)^/^ is sensitive to them all. Thus, det will be needed 
for that theory (see Chapter [5]) . 

There will also be a left inner product and, correspondingly, two sets of MOPRL and MOPUC. 
We discuss this further in Sections 12.11 and 13. li 

Occasionally, for C' vector- valued functions / and g, we will want to consider the scalar 

j 7U^Mkjix)gjix)d^itrix), (1.34) 

which we will denote 

^ d{f{x),^,(x)g(x))ci. (1.35) 



We next turn to module Fourier expansions. A set {v^jjjLi ^ may be infinite) is called 
orthonormal if and only if 

{{^j,^k))R = Sjk'i- (1-36) 

This natural terminology is an abuse of notation since (jl.36p implies orthogonality in {■,-)r but 
not normalization, and is much stronger than orthogonality in (•, ■)r. 

Suppose for a moment that N < oo. For any oi, . . . , ajv G Mi, we can form Ylf=i ^j^j ^'^d, by 
the right multiplication relations ()1.28p . (jl.29p . and (|1.36p . we have 



N N N 



i=i i=i ^ i=i 

We will denote the set of all such Yl!j=i VjO-j by "^(^5^) — it is a vector subspace of Ti of dimension 
(over C) 
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Define for f GH, 

N 

^{^,)U) = Y.^d^3J))R- (1-38) 
i=i 

It is easy to see it is the orthogonal projection in the scalar inner product (•, from Ti to 'H(^^.y 
By the standard Hilbert space calculation, taking care to only multiply on the right, one finds 
the Pythagorean theorem, 

TV 

((/, = ((/ - ^(^,)/, / - ^i^,)f))R + Y,ii^jJ))Uvj^ f))R- (1-39) 

i=i 

As usual, this proves for infinite that 

N 

i=i 

and the convergence of 

N 

J2^j((^^^f))R^''i^Mf) (1-41) 

i=i 

allowing the definition of vr(^^,) and of = Ran7r(^^,) for N = oo. 

An orthonormal set is cahed complete if 'H(ip-) = Ti.. In that case, equality holds in (jl.40p and 
^(^j) (/)=/• 

For orthonormal bases, we have the Parseval relation from p.39p 

oo 

{{fJ))R = T.{{^jJ))Ui^jJ))R (1-42) 
i=i 

and 

oo 

\\f\\l = Y.^^i{{v^,nUiv,j))R)- (1-43) 

1.3 Matrix Mobius Transformations 

Without an understanding of matrix Mobius transformations, the form of the MOPUC Geronimus 
theorem we will prove in Section [3.101 will seem strange-looking. To set the stage, recall that scalar 
fractional linear transformations (FLT) are associated to matrices T = ^) with detT 7^ via 

Without loss, one can restrict to 

det(r) = 1. (1.45) 

Indeed, T fx is a 2 to 1 map of §L(2, C) to maps of C U {00} to itself. One advantage of the 
matrix formalism is that the map is a matrix homomorphism, that is, 



Jtos = It ° fs, 



(1.46) 
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which shows that the group of FLTs is SL(2,C)/{1, —1}. 

While (jl.46p can be checked by direct calculation, a more instructive way is to look at the 
complex projective line, n, f G \ {0} are called equivalent if there is A G C \ {0} so that u = Xv. 
Let [•] denote equivalence classes. Except for [(g)], every equivalence class contains exactly one 
point of the form (^) with z G C. If [(q)] is associated with oo, the set of equivalence classes is 
naturally associated with C U {oo}. /t then obeys 



(1.47) 



KO] 













from which (I1.46P is immediate. 

By Mobius transformations we will mean those FLTs that map B onto itself. Let 



J 



(L48) 



Then [u] = [[fj] with \z\ = 1 (resp. \z\ < 1) if and only if {u, Ju) = (resp. {u,Ju) < 0). From 
this, it is not hard to show that if det(T) = 1, then fx maps D invertibly onto B if and only if 



T^JT = J. 



If T has the form ( " ^ ) , this is equivalent to 



\b\'-\d\' 



ab — cd = 0. 



(1.49) 
(1.50) 



The set of T's obeying det(T) = 1 and ()1.49p is called SU(1, 1). It is studied extensively in [168^ 
Sect. 10.4]. 

The self-adjoint elements of Sl[J(l, 1) are parametrized by a G B via p = (1 — jap)^/^. 



associated to 

Notice that 
and that 



1 A a 
p\a 1 

z + a 



1 + az 



(1.51) 

(1.52) 
(1.53) 



\/z G 



z, 



3 ! a such that Tq,(0) 
namely, a = z. 

It is a basic theorem that every holomorphic bijection of D to D is an fx for some T in §U(1, 1) 
(unique up to ±1). 

With this in place, we can turn to the matrix case. Let A4i be the space of / x / complex 
matrices with the Euclidean norm induced by the vector norm (• 



Let 



{AG Mr. \\A\\ < 1}. 



(1.54) 



We are interested in holomorphic bijections of B; to itself, especially via a suitable notion of FLT. 
There is a huge (and diffuse) literature on the subject, starting with its use in analytic number 
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theory. It has also been studied in connection with electrical engineering filters and indefinite 
matrix Hilbert spaces. Among the huge literature, we mention [H [3l [78l [991 11141 1166| . Especially 
relevant to MOPUC is the book of Bakonyi-Constantinescu [6]. 

Consider Mi © Mi = Mi[2] as a right module over Mi- The A^;-projective line is defined by 
saying [y ] ~ [y']' both in A4/[2] \ {0}, if and only if there exists A E Mi, A invertible so that 



X = X'A, Y = Y'A. 



Let T be a map of A^i[2] of the form 



acting on 7W;[2] by 



T 



A B 
C D 



'X 




'AX + BY' 


Y 




CX + DY 



(1.55) 
(1.56) 
(1.57) 



Because this acts on the left and A equivalence on the right, T maps equivalence classes to them- 
selves. In particular, if CX + D is invertible, T maps the equivalence class of [ ^ ] to the equivalence 



class of 



MX] 
1 



where 



friX] = {AX + B){CX + D)-\ 



(1.58) 



(1.59) 



'X' 


t 


'X' 




J 




1 


1 



So long as CX + D remains invertible, ()1.46p remains true. Let J be the 21 x 21 matrix in / x / 
block form 

Note that (with [f ]^ = [Xtl]) 

it r 

< O'^X^'X <1<^ \\X\\ < 1. 

Therefore, if we define SU(/, /) to be those T's with detT = 1 and 

tUt = J, 

then 

T G SU(/,/) frpi] = as a bijection. 
If T has the form (ll.56p . then (ll.6ip is equivalent to 

A^A - C^C = D^D - B^B = 1, 
A^B = C^D 

(the fourth relation B'^ A = D^C is equivalent to ()1.64p ). 
This depends on 



(1.60) 

(1.61) 
(1.62) 

(1.63) 
(1.64) 



Proposition 1.1. If T = {c d) obeys (ll.6ip and \\X\\ < I, then CX + D is invertible. 
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Proof. ()1.6ip implies that 

T-^ = JTU (1.65) 

-i ~D* )- 

Clearly, ([LHT]) also implies T'^ G §U(^,/). Thus, by ([LM]) for p-^, 

DD^ - CC^ = 1. (1.67) 

This implies first that DD'^ > 1, so D is invertible, and second that 

\\D-'^C\\ < 1. (1.68) 

Thus, ||X|| < 1 implies \\D~'^CX\\ < 1 so l + D-'^CX is invertible, and thus so is D{1 + D-'^CX). 

□ 

It is a basic result of Cartan [IS] (see Helgason |114j and the discussion therein) that 
Theorem 1.2. A holomorphic bijection, g, of Hi to itself is either of the form 

g{X) = friX) (1.69) 

for some P G SU(/,/) or 

g{X) = friX'). (1.70) 

Given a G Aii with ||a|| < 1, define 

/ = (1 - ata)V2, pR = (1 _ «at)i/2. (1.71) 

Lemma 1.3. We have 

ap^ = p^a, a^p^ = p^a\ (1.72) 

Proof. Let / be analytic in D with f{z) = X^^q Taylor series at 2: = 0. Since ||a^a|| < 1, 
we have 

00 

f{a^a) = Y,Cn{a^<^T (1-74) 

n=0 

norm convergent, so a{a^a)^ = {aa^)^a implies 

a/(a^Q) = f{aa^)a, (1-75) 

which implies the first halves of (jl.72p and ()1.73p . The other halves follow by taking adjoints. □ 

Theorem 1.4. Phere is a one-one correspondence between a's in Mi obeying \\a\\ < 1 and positive 
self-adjoint elements o/SU(/,/) via 



p - i ^p""^'' ip^'r'^] n 76) 
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Proof. A straightforward calculation using Lemma [1 . 31 proves that is self-adjoint and TaJTa = J. 

Conversely, if T is self-adjoint, T = and in SU(/,/), then = T ^ A'f = A, = C , so 
(I1.63P becomes 

AA^ - BB^ = 1, (1.77) 

so if 

a = A-^B, (1.78) 

then ()1.77p becomes 

A-\A-^)^ + aa^ = 1. (1.79) 

Since T > 0, A > so pTTO]) implies A = {p^y^, and then (fTTS]) implies B = {p^^^a. 
By Lemma 11.31 

C = B^ = a\p^)-^ = {p^y^a^ (1.80) 
and then (by D = D\ = B, and (fTMIl ) DD^ - CC^ = 1 plus D>0 implies D = (/)"^ □ 
Corollary 1.5. For each a G B;, the map 

ItAX) = (p^)-^(X + a)(l + atx)-i(/) (1.81) 
takes B/ to B;. Its inverse is given by 

/£\X) = fT_jX) = {p^r'iX - a){l - a^X)-\p^). (1.82) 
There is an alternate form for the right side of (jl.Sip . 
Proposition 1.6. The following identity holds true for any X , \\X\\ < 1: 

p^'il + Xa^)-\X + a){p'^y^ = {p'')-\X + a)(l + X)-^ p"^ . (1.83) 
Proof. By the definition of p^ and p^ , we have 

X{p^)-^{1 - a^a) = (p^)-2(l - aa^)X. 
Expanding, using (jl.73p and rearranging, we get 

X(p-^)-2 + a(/)-2atx = [p^y^X + Xa^ip^y^a. 
Adding a{p^y'^ + X{p^)^'^a^ X to both sides and using ()1.73p again, we obtain 

X{p^y^ + a{p^y^ + X{p^y^a^X + a{p^y^a^X 

= {p^y^X + {p^y^a + Xa^p^y^X + Xa^p^y^a, 

which is the same as 

{X + a)(/)-'(l + a^X) = (1 + Xat)(p«)-2(X + a). 
Multiplying by (1 + Xa^y^ and (1 + a^Xy^, we get 

(1 + Xa^y\X + a)(/)-2 = {pRy\X + a)(l + Xy^ 
and the statement follows. □ 
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1.4 Applications and Examples 

There are a number of simple examples which show that beyond their intrinsic mathematical 
interest, MOPRL and MOPUC have wide application. 

(a) Jacobi matrices on a strip 

Let A C be a subset (perhaps infinite) of the z^-dimensional lattice Z'^ and let £^(A) be square 
summable sequences indexed by A. Suppose a real symmetric matrix aij is given for all i,j G A 
with aij = unless |i — j| = 1 (nearest neighbors). Let (3i be a real sequence indexed by i G A. 
Suppose 

suplojjl +sup|/3j| < oo. (1-84) 

i,j i 

Define a bounded operator, J, on ^^(A) by 

{Ju)i = ^ atjUj + PiUi. (1.85) 

j 

The sum is finite with at most 2^ elements. 

The special case A = {1,2, ... } with bi = (5i, ai = ai,j+i > corresponds precisely to classical 
semi-infinite tridiagonal Jacobi matrices. 

Now consider the situation where A' C Ij^^"^ is a finite set with / elements and 

A = {j G : Ji G {1, 2, . . . }; (is, • • • J.) G A'}, (1.86) 
a "strip" with cross-section A'. J then has a block / x / matrix Jacobi form where (7, 5 G A') 

(i?,)75 = %,7)' (7 = ^): (1-87) 

= a(i,7)(i,<5), (t/*^), (1-88) 

{Ai)^s = a(i,7)(i+i,5)- (1.89) 

The nearest neighbor condition says {Ai)^^ = if 7 7^ (5. If 

> (1.90) 

for all i, 7, then Ai is invertible and we have a block Jacobi matrix of the kind described in Section [212] 
below. 

By allowing general Ai,Bi, we obtain an obvious generalization of this model — an interpretation 
of general MOPRL. 

Schrodinger operators on strips have been studied in part as approximations to Z'^; see |3 H 195 1 
11301 1134] I15H[T6^ . Prom this point of view, it is also natural to allow periodic boundary conditions 
in the vertical directions. Furthermore, there is closely related work on Schrodinger (and other) 
operators with matrix- valued potentials; see, for example, PlHESlEelETllSllEQllMl 1971 [165]. 
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(b) Two-sided Jacobi matrices 

This example goes back at least to Nikishin [153]. Consider the case v = 2, A! = {0, 1} C Z, and 
A as above. Suppose (jl.90p holds, and in addition, 

a{i,o){i,i) > 0, (1.91) 
a(i,o){j,i) = 0, i = 2,3, .... (1.92) 

Then there are no links between the rungs of the "ladder," {1, 2, . . . } x {0, 1} except at the end 
and the ladder can be unfolded to Z! Thus, a two-sided Jacobi matrix can be viewed as a special 
kind of one-sided 2x2 matrix Jacobi operator. 

It is known that for two-sided Jacobi matrices, the spectral theory is determined by the 2x2 
matrix 



11 



where d^jiki is the measure with 



{6u,{J-\)-Hi) = j (1.94) 

but also that it is very difficult to specify exactly which d^i correspond to two-sided Jacobi matrices. 

This difficulty is illuminated by the theory of MOPRL. By Favard's theorem (see Theorem l2.1ip . 
every such d/i (given by (11.93P and positive definite and non-trivial in a sense we will describe in 
Lemma [2T]) yields a unique block Jacobi matrix with Aj > (positive definite). This d/j, comes 
from a two-sided Jacobi matrix if and only if 

(a) Bj is diagonal for j = 2,3, . . . . 

(b) Aj is diagonal for j = 1, 2, . . . . 

(c) Bj has strictly positive off-diagonal elements. 
These are very complicated indirect conditions on dfil 



(c) Banded matrices 

Classical Jacobi matrices are semi-infinite symmetric tridiagonal matrices, that is, 

Jkm = if|fc-m|>l (1.95) 

with 

Jkm>0 if|/c-m| = l. (1.96) 
A natural generalization are {21 + l)-diagonal symmetric matrices, that is, 

Jkni = if\k-m\>l, (1.97) 
Jfc„>0 if|A;-m|=/. (1.98) 

Such a matrix can be partitioned into / x I blocks, which is tridiagonal in block. The conditions 
(jl.97p and (|1.98p are equivalent to Ak E C, the set of lower triangular matrices; and conversely, 
j4fc E £, with A}^,Bk real (and 5^ symmetric) correspond precisely to such banded matrices. This 
is why we introduce type 3 MOPRL. 

Banded matrices correspond to certain higher-order difference equations. Unlike the second- 
order equation (which leads to tridiagonal matrices) where every equation with positive coefficients 
is equivalent via a variable rescaling to a symmetric matrix, only certain higher-order difference 
equations correspond to symmetric block Jacobi matrices. 
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(d) Magic formula 

In [53], Damanik, Killip, and Simon studied perturbations of Jacobi and CMV matrices with 
periodic Jacobi parameters (or Verblunsky coefficients). They proved that if A is the discriminant 
of a two-sided periodic Jq, then a bounded two-sided J has A(J) = 5^ -|- 5"^ {{Su)n = Un+i) if 
and only if J lies in the isospectral torus of Jq. They call this the magic formula. 

This allows the study of perturbations of the isospectral torus by studying A(J) which is a 
polynomial in J of degree p, and so a 2p -|- 1 banded matrix. Thus, the study of perturbations of 
periodic problems is connected to perturbations of + S^^ as block Jacobi matrices. Indeed, it 
was this connection that stimulated our interest in MOPRL, and [33j uses some of our results here. 

(e) Vector-valued prediction theory 

As noted in Section 11.11 both prediction theory and filtering theory use OPUC and have natural 
MOPUC settings that motivated much of the MOPUC literature. 

2 Matrix Orthogonal Polynomials on the Real Line 
2.1 Preliminaries 

OPRL are the most basic and developed of orthogonal polynomials, and so this chapter on the 
matrix analogue is the most important of this survey. We present the basic formulas, assuming 
enough familiarity with the scalar case (see [231 EH 11671 11761 11841 1185j ) that we do not need to 
explain why the objects we define are important. 

2.1.1 Polynomials, Inner Products, Norms 

Let dfj, he an I X I matrix-valued Hermitian positive semi-definite finite measure on M normalized 
by /x(M) = 1 G Ml- We assume for simplicity that /x has a compact support. However, many of 
the results below do not need the latter restriction and in fact can be found in the literature for 
matrix-valued measures with unbounded support. 
Define (as in ([120])) 

{{f,g))R = I f{x)Ufi{x)g{x), ||/||h = (TY((/,/))^)1/2, f,geV, 

{{f,g))L = l9{x)df,{x)f{x)l ||/||z. = (Tr((/,/))^)V2, 

Clearly, we have 

((/,9))L = ((5,/))i?o {{f,9))i = {{9J))L, (2.1) 

{{f,9))L = {{9^,P))R, II/IIl = ||/^lb. (2.2) 

As noted in Section 11.21 we have the left and right analogues of the Cauchy inequality 
|Tr((/,5))ij| < UWrMr, |Tr((/,g))L| < ||/||l Ml- 

Thus, ||-||_R and ||-||l are semi-norms in V. Indeed, as noted in Section [L2l they are associated 
to an inner product. The sets {/ : ||/||/j = 0} and {/ : ||/||l = 0} are linear subspaces. Let Vji 
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be the completion of V /{f : ||/||_r = 0} (viewed as a right module over Mi) with respect to the 
norm ||-||_r. Similarly, let Vl be the completion oiV /{f : = 0} (viewed as a left module) with 

respect to the norm 

The set V defined in Section 11.21 is a linear space. Let us introduce a semi-norm in V by 



Let Vo C V be the linear subspace of all polynomials such that |/| = and let Vqo be the completion 
of the quotient space V/Vq with respect to the norm | • |. 

Lemma 2.1. The following are equivalent: 

(1) ll/lliJ > for every non-zero f GV. 

(2) For all n, the dimension in Vr of the set of all polynomials of degree at most n is {n + 1)P . 

(3) II/IIl > for every non-zero / G "P. 

(4) For all n, the dimension in Vl of the set of all polynomials of degree at most n is {n + 1)1'^. 

(5) For every non-zero v G V, we have that \v\ ^ 0. 

(6) For all n, the dimension in Voo of all vector-valued polynomials of degree at most n is {n + \)l. 
The measure dfi is called non-trivial if these equivalent conditions hold. 

Remark. If / = 1, these are equivalent to the usual non- triviality condition, that is, supp(^) is 
infinite. For / > 1, we cannot define triviality in this simple way, as can be seen by looking at the 
direct sum of a trivial and non-trivial measure. In that case, the measure is not non-trivial in the 
above sense but its support is infinite. 

Proof. The equivalences (1) 44> (2), (3) <^ (4), and (5) <^ (6) are immediate. The equivalence (1) 
O (3) follows from ()2.2p . Let us prove the equivalence (1) 44> (5). Assume that (1) holds and let 
V G V he non-zero. Let / € A4i denote the matrix that has v as its leftmost column and that has 
zero columns otherwise. Then, ^ ||/||r = Tr((/7 = 1^1^ and hence (5) holds. Now assume 
that (1) fails and let / G "P be non-zero with ||/||_r = 0. Then, at least one of the column vectors 
of / is non-zero. Suppose for simplicity that this is the first column and denote this column vector 
by V. Let t G A4i be the matrix tij = 5ii5ji; then we have 



Throughout the rest of this chapter, we assume the measure dfi to be non-trivial. 

2.1.2 Monic Orthogonal Polynomials 

Lemma 2.2. Let dfj, be a non-trivial measure. 
(i) There exists a unique monic polynomial P^ of degree n, which minimizes the norm \\P^\\ji. 




(2.3) 



WfU = {{fJ))R = 0^0 = Tv{t*{{f,f))Rt) = \v\' 



and hence (5) fails. 



□ 
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(iv) The polynomial can he equivalently defined as the monic polynomial of degree n which 
satisfies 

{{P^^,f))L = foranyfeV, deg/ < n. (2.5) 

(v) One has P^ix) = Pn{x)'^ for all x and 

{{P^,P^))r = {{P^,P^))l. (2.6) 

Proof. As noted, V has an inner product (•, so there is an orthogonal projection vr^^^ onto Vn 
discussed in Section 11.21 Then 

Pi^(x) = x'^-vrS(x"). (2.7) 

As usual, in inner product spaces, this uniquely minimizes x^ — Q over all Q G Vn-i- It clearly 
obeys 

Tr(((P„^/))^) = (2.8) 
for all / G Vn-i- But then for any matrix a, 

TTi{{P^J))Ra) = Tv{{{P^Ja))R) = 



so ([23]) holds. 

This proves (i) and (ii). (iii) and (iv) are similar. To prove (v), note that P^(x) = P^(x)t 
follows from the criteria (fTi]) . (p3|) . The identity ^M) follows from (fOj) . □ 



Lemma 2.3. Let he non-trivial. For any monic polynomial P, we have det{{P,P))R / and 
det((P,P))L /O. 

Proof. Let P be a monic polynomial of degree n such that ((P, P))r has a non-trivial kernel. Then 
one can find a G A4i, a ^ 0, such that a^{{P,P))jia = 0. It follows that ||Pa||/j = 0. But since 
P is monic, the leading coefficient of Pa is a, so Pa ^ 0, which contradicts the non-triviality 
assumption. A similar argument works for ((P, P))/,. □ 

By the orthogonality of Qn — Pn to P^ for any monic polynomial Qn of degree n, we have 

((Qn, Qn))R = HQ - P^, Q - P^))r + {{P^, P^))r (2.9) 

and, in particular, 

{{P^,P^))R<{{Qn,Qn))R (2.10) 

with (by non-triviality) equality if and only if Qn = P^- Since Tr and det are strictly monotone 
on strictly positive matrices, we have the following variational principles ( ()2.1ip restates (i) of 
Lemma 12.211 : 



Theorem 2.4. For any monic Qn of degree n, we have 

\\Qn\\R>\\P^\\R, (2.11) 

det((Q„, Qn))R > det((P„^, P^))r (2.12) 
with equality if and only if Pn = Qn ■ 
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2.1.3 Expansion 

Theorem 2.5. Let dfi be non-trivial. 

(i) We have 

{{Pi',P^))R = JnSkn (2.13) 

for some positive invertihle matrices 7„. 

(ii) {-P^I^^Q are a right-module basis forVni indeed, any f £ Vn has a unique expansion, 

n 

f = Y.Ppfj'- (2.14) 

i=o 

Indeed, essentially by (jl.38p . 

fp = ^7'{{Pj'j))R- (2-15) 
Remark. There are similar formulas for {{■,■)) l- By (j2.6p . 

{{Pt,P^))L=ln6kn (2.16) 

(same jn, which is why we use and not 7^). 

Proof, (i) (I2.13P for n < A; is immediate from (I2.5p and for n > A; by symmetry. 7n > follows 
from (jl.22p . By Lemma 12.31 det(7„) 7^ 0, so 7„ is invertible. 

(ii) Map {Mir+^ to Vn by 

n 

{ao, . . . , a„) 1-^ ^ Pf-Uj = X{ao, ««). 
3=0 



By mM, 

"'■3 - ij j 

so that map is one-one. By dimension counting, it is onto. □ 



«i = 7,- H{Pf^ ^("0, . . . , an))) 



2.1.4 Recurrence Relations for Monic Orthogonal Polynomials 

Denote by Cn (resp. ) the coefficient of in P^ix) (resp. P^(x)), that is, 

Pn{x) = x^l + CnX""'^ + lower order terms, 
P^(x) = x^l + CnX"^'^ + lower order terms. 

Since P^(x)^ = P^(x), we have (Cn)^ = Cn - Using the parameters 7^ of ()2.13p and j Cn can 
write down recurrence relations for Pn{x), P^(x). 

Lemma 2.6. (i) We have a commutation relation 

7n-l(e - Ctl) = iCf: - Ctlhn-l. (2.17) 

(ii) We have the recurrence relations 

xP^ix) = P„^+i(x) + P^{x){C^ - C^+i) + P„_i(x)7-ii7n, (2.18) 

xP^ix) = PiVl(^) + (Cn - Cn+l)^n (^) +7n7nil^n-l(^)- (2-19) 
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Proof, (i) We have 

- xP^_i{x) = {Cn - Cn-i)^;""^ + lower order terms 

and so 

— {{Pn ~ ^Pn-1^ Pn-l))R 

= {{P^,P^_,))r - {{xPti,Pti))R 
= -{{xPti,Ptl))R 
= -{{Pti,xPti))R 

— {{Pn~l^ Pn ~ xPn-lflR 

= {{Pti.x--\C^-Cti)))R 

= {{Ptl.X^-^))R{C^-Ctl) 
= 7n-li(n - Cn-l)- 

(ii) By Theorem 12. 51 

xP^ix) = P^^^{x)Cn+l + P^{x)Cn + P^_^{x)Cn-l + ■■■+ P^Cq 

with some matrices Co, ... , Cn+i- It is straightforward that C„+i = 1 and Cn = Cn ~ Cn+i- ^he 
orthogonahty property (j2.4p . we find Co = • • • = C„_2 = 0- Finahy, it is easy to calculate C„_i: 

7„ = {{P^,xPti))R = {{xP^,Pti))R 

= {{Pn+l + PniCn ~ Cn+l) + Pn-l^ n-l , Pn-l)) R 
= Ci_i7n-1 

and so, taking adjoints and using self-adjointness of 7j, C„_i = 7^j';i7n. This proves (j2.18p : the 
other relation (j2.19p is obtained by conjugation. □ 

2.1.5 Normalized Orthogonal Polynomials 

We call £ V a right orthonormal polynomial if degp^ < n and 

{{Pn, f))R = for every feV with deg f < n, (2.20) 
{{p^,P^))r = 1. (2.21) 

Similarly, we call & V a left orthonormal polynomial if degp^ < n and 

{{Pn, f))L = for every / G P with deg / < n, (2.22) 
((Pn,Pn))L = l- (2.23) 
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Lemma 2.7. Any orthonormal polynomial has the form 

p^ix) = P^ix){{P^, P^Tn'^an, Pni^) = TniiP^, P^))l'/^ P^ {x) (2.24) 

where an,Tn E ^Al are unitaries. In particular, degp^ = degp^ = n. 

Proof. Let Kn be the coefficient of in p^. Consider the polynomial q[x) = P^{x)Kn — p^{x), 
where is the monic orthogonal polynomial from Lemma 12.21 Then degq' < n and so from l\2.4\\ 
and (j2.20p . it follows that {{q,q))R = and so q{x) vanishes identically. Thus, we have 

1 = {{p^,P^))r = KliiP^, P^llRKn (2.25) 
't\-li<'-l _ II jdR d-Rw „ „y^A li^ Zi't _ // D-R i3-R\\-l 



and so det(i^„) / 0. From we get (K^)-!^-! ^ ^^pK^ pHj^^^ ^nd so KnKl = {{P^, P^))^' 

From here we get Kn = {{P^ , Pn)) r^"^ '^n with a unitary (T„. The proof for p^ is similar. □ 

By Theorem \2.b\ the polynomials p^ form a right orthonormal module basis in Vr. Thus, for 
any / E Vr, we have 

oo 

f{x) = = '^^^^ (2.26) 

m=0 

and the Parseval identity 

oo 

Y,^rUmfl) = \\f\?R (2.27) 

m=0 

holds true. Obviously, since / is a polynomial, there are only finitely many non-zero terms in (j2.26p 
and (12:271) . 

2.2 Block Jacobi Matrices 

The study of block Jacobi matrices goes back at least to Krein [133j . 
2.2.1 Block Jacobi Matrices as Matrix Representations 

Suppose that a sequence of unitary matrices 1 = uo, cJi, ct2, . . . is fixed, and p^ are defined according 
to (j2.24p . As noted above, p^ form a right orthonormal basis in Vr. 

The map /(x) i— > xf{x) can be considered as a right homomorphism in Vr. Consider the matrix 
Jnm of this homomorphism with respect to the basis p^, that is, 

Jnm = {{p^-l,Xp^_^))R. (2.28) 

Following Killip-Simon |128j and Simon |167|, 11681 I176j . our Jacobi matrices are indexed with 
n = 1,2,... but, of course, p„ has n = 0, 1, 2, . . . . That is why (j2.28p has n — 1 and m — 1. 

As in the scalar case, using the orthogonality properties ofp^, we get that Jnm = if |n— m| > 1. 
Denote 

Bn = Jnn = {{Pn-1, XPn-l)) R 

and 

An = Jn,n+l = <^I+l,n = {{Pn-1^ ^Pn)) R- 
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Then we have 



J 



/Si ^1 •••\ 
4 B2 A2 •• 



Al Bs 



V 



(2.29) 



Applying (j2.26p to f{x) = xp^{x), we get the recurrence relation 

xPnix) = p^+iix)Al^-^ + p^{x)Bn+i +p^„i(x)>l„, n = 1, 2, , 



(2.30) 



If we set p^i{x) = and = Ij the relation (j2.30p also holds for n = 0. By (j2.2p . we can always 
pick so that for x real, p^i^) = Pni^V ^ thus for complex z, 



p^iz)=p^{z)^ 

by analytic continuation. By conjugating (j2.30|) . we get 

XPnix) = An+ip^_^^{x) + Bn+lPn{x) + A]jP.^_i(x), n = 0, 1,2, . . . . 

Comparing this with the recurrence relations (j2.18p . (j2.19p . we get 

, t -1/2 1/2 r, t 1/2 /aJ? a_R\ -1/2 

= <-l7n-l 7n f^n, = <_i7„_l(Cn-l " C„ )7„-l C^n-l- 



(2.31) 



(2.32) 



(2.33) 



In particular, det An 7^ for all n. 

1/2 1/2 

Notice that since (T„ is unitary, |det((Tn)| = 1, so (|2.33p implies det(7n ) = det(7^_]^)|det(74„)| 
which, by induction, implies that 



det((P„«,Pi^)) = |det(Ai---A„ 



(2.34) 



Any block matrix of the form (j2.29p with i?„ = B^ and det An 7^ for all n will be called a 
block Jacobi matrix corresponding to the Jacobi parameters An and 

2.2.2 Basic Properties of Block Jacobi Matrices 

Suppose we are given a block Jacobi matrix J corresponding to Jacobi parameters An and Bn, 
where Bn = Bn and det An for each n. 

Consider the Hilbert space T-Cv = ^^(Z+, C') (here Z+ = {1, 2, 3, . . . }) with inner product 

00 

n=l 

and orthonormal basis {ek,j}k&+,i<j<i, where 

and {vj}i<j<i is the standard basis of C'. J acts on 7iv via 

iJf)n = Al^Jn^l + Bnfn + Anfn+l, f & (2.35) 
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(with /o = 0) and defines a symmetric operator on this space. Note that using invertibihty of the 
A„'s, induction shows 

span{efcj : I < k < n, I < j < 1} = span{ j'^^^eij : I < k < n, I < j <l} (2.36) 

for every n > 1. We want to emphasize that elements of Ti^ and 7i are vector- valued and matrix- 
valued, respectively. For this reason, we will be interested in both matrix- and vector-valued 
solutions of the basic difference equations. 

We will consider only bounded block Jacobi matrices, that is, those corresponding to Jacobi 
parameters satisfying 

sup TV(4A„ + BlBn) < oo. (2.37) 

n 

Equivalently, 

sup(P„|| + ||5„||) < oo. (2.38) 

n 

In this case, J is a bounded self-adjoint operator. This is equivalent to fi having compact support. 

We call two Jacobi matrices J and J equivalent if there exists a sequence of unitaries Un £ Mi, 
n > 1, with ui = 1 such that Jnm = UnJnmUm- From Lemma 12.71 it is clear that if p^, 
are two sequences of normalized orthogonal polynomials, corresponding to the same measure (but 
having different normalization), then the Jacobi matrices Jnm = {{Pn-n xp^_i)) r and Jnm — 
{{Pn-i-,xp^_i))R are equivalent (n„ = a^_i(T„_i). Thus, 

Bn = ulBnUn, i„ = n]j^„n„+i . (2.39) 

Therefore, we have a map 

$ : ^ ^{J : Jmn = {{Pn-i,xp^-i))R, Pn ^2 40) 

correspond to dfx for some normalization} 

from the set of all Hermitian positive semi-definite non-trivial compactly supported measures to 
the set of all equivalence classes of bounded block Jacobi matrices. Below, we will see how to invert 
this map. 



2.2.3 Special Representatives of the Equivalence Classes 

Let J be a block Jacobi matrix with the Jacobi parameters An, Bn- We say that J is: 

• of type 1, if An > for all n; 

• of type 2, if A1A2 • • • > for all n; 

• of type 3, if An € C for all n. 

Here, C is the class of all lower triangular matrices with strictly positive elements on the diagonal. 
Type 3 is of interest because they correspond precisely to bounded Hermitian matrices with 2/ + 1 
non- vanishing diagonals with the extreme diagonals strictly positive; see Section [1.4( c). Type 2 is 
the case where the leading coefficients of p^ are strictly positive definite. 

Theorem 2.8. (i) Each equivalence class of block Jacobi matrices contains exactly one element 
each of type 1, type 2, or type 3. 
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(ii) Let J be a block Jacobi matrix corresponding to a sequence of polynomials p^ as in (j2.24p . 
Then J is of type 2 if and only if an = 1 for all n. 

Proof. The proof is based on the following two well-known facts: 

(a) For any t G Mi with det(t) 7^ 0, there exists a unique unitary u E Mi such that tu is Hermitian 
positive semi-definite: tu > 0. 

(b) For any t G Mi with det(t) 7^ 0, there exists a unique unitary u G Mi such that tu £ C. 

We first prove that every equivalence class of block Jacobi matrices contains at least one element 
of type 1. For a given sequence An, let us construct a sequence tii = 1, ^2, M3, . . . of unitaries such 
that UnAnUn+i > 0. By the existence part of (a), we find U2 such that A1U2 > 0, then find M3 
such that U2A2U3 > 0, etc. This, together with (12.391) . proves the statement. In order to prove the 
uniqueness part, suppose we have An > and UnAnUn+i > for all n. Then, by the uniqueness 
part of (a), Ai > and A1U2 > imply U2 = 1; next, > and u\A2Us = A2U3 > imply 
U3 = 1, etc. 

The statement (i) concerning type 3 can be proven in the same way, using (b) instead of (a). 

The statement (i) concerning type 2 can be proven similarly. Existence: find U2 such that 
A1U2 > 0, then ii3 such that {AiU2){u\A2U3) = A1A2U3 > 0, etc. Uniqueness: if Ai . . .An > 
and Ai - ■ ■ AnUn+i > 0, then Un+i = 1. 

By ()2.33p . we have A1A2 ■ ■ ■ An = 7n o"n and the statement (ii) follows from the positivity of 

In- □ 

We say that a block Jacobi matrix J belongs to the Nevai class if 

Bn — > and Al^An ^ 1 as n — > 00. 

It is clear that J is in the Nevai class if and only if all equivalent Jacobi matrices belong to the 
Nevai class. 

Theorem 2.9. // J belongs to the Nevai class and is of type 1 or type 3, then An ^1 as n ^ 00. 

Proof. If J is of type 1, then An An = — > 1 clearly implies A„ — > 1 since square root is continuous 
on positive Hermitian matrices. 

Suppose J is of type 3. We shall prove that An ^ 1 by considering the rows of the matrix An 
one by one, starting from the /th row. Denote {An)jk = o^^j^k- have 

{AiAnh = {a[f)^ ^ 1, and so a[f ^ 1. 

Then, for any k < I, we have 

{AlAn)ik = a^fa^"}! ^ 0, and so a\'^^ 0. 
Next, consider the (/ — l)st row. We have 

(4A„)z-i,/-i = (a;Vi)' + l«&-iP-l 

(n) 

and so, using the previous step, a^ -^ ^ ^ 1 as n — > 00. Then for all A; < Z — 1, we have 

[Ai,An)i-i,k - az-i,fc + az,/-i ^k ^ ^ 

(n) 

and so, using the previous steps, az-i,fc 0. Continuing this way, we get — > 5j^k as required. □ 
It is an interesting open question if this result also applies to the type 2 case. 
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2.2.4 Favard's Theorem 

Here we construct an inverse of the mapping $ (defined by (|2.4U|) ). Thus, ^ sets up a bijection 
between non-trivial measures of compact support and equivalence classes of bounded block Jacobi 
matrices. 

Before proceeding to do this, let us prove: 
Lemma 2.10. The mapping $ is injective. 

Proof. Let // and jl be two Hermitian positive semi-definite non-trivial compactly supported mea- 
sures. Suppose that = 

Let and p^ be normalized orthogonal polynomials corresponding to /x and /i. Suppose that 
the normalization both for and for p^ has been chosen such that fj^ = 1 (see (|2.24p ). that 
is, type 2. From Lemma 12.81 and the assumption $(//) = $(/i) it follows that the corresponding 
Jacobi matrices coincide, that is, {{PnT^Pm))R — {{Pnj^PnJjR for all n and m. Together with the 
recurrence relation (|2.30p this yields p^ = p^ for all n. 

For any n > 0, we can represent x" as 



k=0 k=0 



The coefficients C^"^ and cj^^ are completely determined by the coefficients of the polynomials p^ 

and p^ and so C^"^ = cj^^ for all n and k. 
For the moments of the measure n, we have 



/ 



x^d^,{x) = {{l,x-))n = = {{1,1)) RCt^ = Ct^. 

k=0 



Since the same calculation is valid for the measure /2, we get 
for all n. It follows that 



f{x)dfi{x)g{x) = J f{x)dfi{x)g{x) 
for all matrix-valued polynomials / and g, and so the measures /i and jl coincide. □ 

We can now construct the inverse of the map ^. Let a block Jacobi matrix J be given. By 
a version of the spectral theorem for self-adjoint operators with finite multiplicity (see, e.g., [21 
Sect. 72]), there exists a matrix- valued measure dfj, with 



(eij,/(J)ei,fe)H„ = J /(x)d^j-fc(x) (2.41) 

and an isometry 

R-.Hv^ L^(R,dfi;d) 
such that (recall that {vj} is the standard basis in C') 

[Reij]{x) = Vj, l<j<l, (2.42) 
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and, for any g £ TC^j, we have 

{RJg){x) =x{Rg){x). (2.43) 

If the Jacobi matrices J and J are equivalent, then we have J = U*JU for some U = ©J^i^n, 
ui = 1. Thus, 

{ei,j,f{J)ei^k)nv = {Ueij,f{J)Uei^k)'Hv = {ei,j, f{J)ei,k)Hv 
and so the measures corresponding to J and J coincide. Thus, we have a map 

{J: J is equivalent to J} fi (2.44) 

from the set of all equivalence classes of bounded block Jacobi matrices to the set of all Hermitian 
positive semi-definite compactly supported measures. 

Theorem 2.11. (i) All measures in the image of the map ^ are non- degenerate. 
(ii) ^ o^i = id. 
(ni) ^ = id. 

Proof, (i) To put things in the right context, we first recall that is a norm (rather than a 

semi-norm), whereas | • | on V (cf. (12. 3p ) is, in general, a semi-norm. Using the assumption that 
det(^fc) 7^ for all k (which is included in our definition of a Jacobi matrix), we will prove that | • | 
is in fact a norm. More precisely, we will prove that \p\ > for any non-zero polynomial p G V; by 
Lemma l2.1l this will imply that fi is non-degenerate. 

Let p € V he a non-zero polynomial, degp = n. Notice that ()2.42p and (12.431) give 

[Rj''eij]{x) =x''vj (2.45) 

for every k > and 1 < j < I. This shows that p can be represented as p = Rg, where g = 
Ylk=o '^^fk, and /o, . . . , /n are vectors in Hv such that (/j, ej^klHv = for alH = 0, . . . , n, j > 2, 
k = 1,...,/ (i.e., the only non-zero components of fj are in the first C' in T~iv). Assumption 
degp = n means /„ / 0. 

Since R is isometric, we have \p\ = WgW-Hv-i and so we have to prove that g ^ 0. Indeed, suppose 
that g = 0. Using the assumption det{Ak) ^ and the tri-diagonal nature of J, we see that 
Y^k=o fk = ^ yields fn = 0, contrary to our assumption. 

(ii) Consider the elements Ren,k £ L^(M, d//; C'). First note that, by (12.36P and (j2.45p . Ren,k is 
a polynomial of degree at most n — 1. Next, by the unitarity of R, we have 

{Ren,k,Rem,j)L2{R,d/_i;Cl) =Sm,nh,j- (2.46) 

Let us construct matrix-valued polynomials qn{x), using Ren,i, Ren,2, ■ ■ ■ , Ren,i as columns of 
Qn-iix): 

[qn-l{x)]j^k = [Ren,k{x)]j. 

We have degg„ < n and {{qm,Qn))R = Sm,n^'i the last relation is just a reformulation of (|2.46p . 
Hence the q^s are right normalized orthogonal polynomials with respect to the measure dfi. We 
find 

Jnm — [(Cnj') <^6m,fc)Wi,]l<j,fc</ 

= [(i?e„j, RJem,k) L^{R^dii;C')]'i-<j,k<l 
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= [{Re.n,j,xRem,k)L'^{R4ix;0)\l<j,k<l 

= [(ten-l(2;)].j,a;[gm-l(a;)].,fc)L2(iR_rf^.C!)]l<j,fc<Z 

= {{qn-l,Xqm-l))R 

as required. 

(iii) Follows from (ii) and from Lemma 12.101 □ 

2.3 The m-Function 

2.3.1 The Definition of the m-Function 

We denote the Borel transform of dji by m: 

m{z) = [ lmz>0. (2.47) 

J X- z 

It is a matrix- valued Herglotz function, that is, it is analytic and obeys lmm{z) > 0. For in- 
formation on matrix-valued Herglotz functions, see [98] and references therein. Extensions to 
operator- valued Herglotz functions can be found in |94j . 

Lemma 2.12. Suppose dfi is given, p^ are right normalized orthogonal polynomials, and J is the 
associated block Jacobi matrix. Then, 



m{z) = {{p^,{x-z)-'p^))R (2.48) 

and 

m{z) = {ei,,{J - z)-'ei,)n^. (2-49) 

Proof. Since Pq = 1, (j2.48p is just a way of rewriting the definition of m. The second identity, 

(j2.49p , is a consequence of (|2.4ip and Theorem I2.11l fiii) . □ 



2.3.2 Coefficient Stripping 

If J is a block Jacobi matrix corresponding to invertible A^s and Hermitian B^s, we denote the 
/c-times stripped block Jacobi matrix, corresponding to {Ak+rn Bk+n}n>i, by J^^\ That IS, 

/Bk+i Ak+i • • A 
4+1 Bk+2 Ak+2 ■■■ 

~ Bk+3 ■■■ ■ 

\ ■ •••/ 

The m-function corresponding to J^^^ will be denoted by m^''\ Note that, in particular, J^^^ = J 
and m^"-* = m. 

Proposition 2.13. Let J be a block Jacobi matrix with (Tess(</) ^ [c^)^]- Then, for every e > 0, 
there is >0 such that for k > ko, we have that a{J^^'^) C [a — e,b + e]. 

Proof. This is an immediate consequence of (the proof of) [42^ Lemma 1]. □ 
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Proposition 2.14 (Due to Aptekarev-Nikishin [3]). We have that 
for Im z > and k > 0. 

Proof. It suffices to liandle tlie case k = 0. Given (j2.49p . tliis is a special case of a general formula 
for 2x2 block operator matrices, due to Schur |163j . that reads 

'A _ f (A-BD-^C)-^ -A-^B{D -CA-^B)-^" 

CD) ~ \-D-^C{A- BD-^C)-^ {D -CA-^B)-^ 

and which is readily verified. Here A = B\ — z, B = Ai, C = A\, and D = J^^-* — z. □ 

2.4 Second Kind Polynomials 

Define the second kind polynomials by q^i{z) = — 1, 

q^{z)=j n = 0,l,2,.... 

Jm z-x 

As is easy to see, for n > 1, (7^ is a polynomial of degree n — 1. For future reference, let us display 
the first several polynomials and q^: 

p^i(x)=0, p^(x) = l, (x) = (x - V> (2.50) 

Q^l(x) = -1, q^{x)=0, qf{x)=A^\ (2.51) 
The polynomials q^ satisfy the equation (same form as (|2.30p ) 

xq^{x)=q^^^{x)Al^^ + q^{x)Bn+i+q^_^{x)An, n = 0,1,2,.... (2.52) 

For n = 0, this can be checked by a direct substitution of (j2.5ip . For n > 1, as in the scalar case, 
this can be checked by taking (|2.30p for x and for z, subtracting, dividing by x — z, integrating 
over d/i, and taking into account the orthogonality relation 



j dij{x)p^{x) =0, n > 1. 



Finally, let us define 

^P^{z)=q^{z)+m{z)p^(z). 
According to the definition of q^, we have 

i.^{z) = {{f,,p^))n, f.{x) = {x--z)-\ 

By the Parseval identity, this shows that for all Imz > 0, the sequence ipni^) is in that is, 

00 

j;Tr(V^(z)V^(z)) <oo. (2.53) 

n=0 

In the same way, we define q^i{z) = —1, 



qi^iz) = I ^M-^dM(x) , n = 0, 1, 2, 



and Tp^{z) = q^{z) +p^{z)m{z). 
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2.5 Solutions to the Difference Equations 

For Im z > 0, consider the solutions to the equations 



ZUn{z) = '^Um{z)Jmn, n = 2, 3, . . . , (2.54) 

m=l 

oo 

ZVniz) = JnmVmiz), n = 2,3, . . . . (2.55) 



m=l 



Clearly, Un{z) solves (j2.54p if and only if Vn{z) = (m„(z))"^ solves (|2.55p . In the above, we normally 
intend z as a fixed parameter but it then can be used as a variable. That is, z is fixed and Un{z) 
is a fixed sequence, not a z-dependent function. A statement like Vn{z) = (ti„(z))^ means if Un is a 
sequence solving (|2.54p for z = zq, then w„ obeys (I2.55P for z = zq. Of course, if Un{z) is a function 
of z in a region, we can apply our estimates to all z in the region. For any solution {uniz)}"^^^ of 
(f23i]) . let us define 

uq{z) = zui{z) — ui{z)Bi — U2{z)A\. (2.56) 

With this definition, the equation ()2.54p for n = 1 is equivalent to uo{z) = 0. In the same way, we 
set 

vo{z) = zvi{z) - Bivi{z) - AiV2{z). 

Lemma 2.15. Let Imz > and suppose {un{z)}'^^Q solves (|2.54p {for n > 2) and (j2.56p and 

belongs to Then 

oo 

(Imz)^TV(u„(z)tu„(z)) = -ImTr(ni(z)no(z)^). (2.57) 

n=l 

In particular, Un{z) = ap^_i{z) is in l"^ only if a = 0. 

Proof Denote s„ = Tr(n„(z)ylJ,_iU„_i(z)t). Here = 1- Multiplying (f23^ forn > 2 and ([236]) 
for n = 1 by Un{z)^ on the right, taking traces, and summing over n, we get 



N N N N 

Z 

71=1 n=l n=l n=l 



^ TT{Uniz)Uniz)'') = ^ + ^ Tr(u„(2;)fi„+iU„(z)^) + ^ Sn . 



Taking imaginary parts and letting — > cxd, we obtain ()2.57p since the middle sum is real and the 
outer sums cancel up to boundary terms. Applying (j2.57p to Un{z) = ap^_i{z), we get zero in the 
right-hand side: 

oo 

(Imz)^Tr(apli(z)p^_i(z)tat) = q 

n=l 

and hence q = since p^ = 1. □ 

Theorem 2.16. Let Imz > 0. 

(i) Any solution {un(z)}^o '^f (!^-54p (for n > 2) can be represented as 

Uniz) = ap^_^{z) + bq^_^{z) (2.58) 
for suitable a,b £ Aii. In fact, a = ui{z) and b = —uo{z). 
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(ii) A sequence (j2.58p satisfies (j2.54p for all n > 1 if and only ifb = 0. 

(iii) A sequence (|2.58|) belongs to i"^ if and only ifun{z) = ciIj^_i{z) for some c G Aii. Equivalently, 
a sequence (12.580 belongs to l"^ if and only if ui{z) -\- UQ{z)m{z) = 0. 

Proof, (i) Let Un{z) be a solution to ()2.54|) . Consider 

Un{z) = Un{z) - Ui{z)p^_i{z) + Uo{z)q^_^{z) . 

Then Un{z) also solves (I2.54p and uq{z) = ui{z) = 0. It follows that Un{z) = for all n. This 
proves (i). 

(ii) A direct substitution of (|2.58p into (I2.54|) for n = 1 yields the statement. 

(iii) We already know that cipn-i is an £^ solution. Suppose that Un{z) is an solution to 
(f23i]l . Rewrite (f238]) as 

Uniz) = {a- bm{z))p^^^{z) + biP^_^{z). 

Since tp^ is in £^ and cp^ is not in we get a = bm{z), which is equivalent to ui{z) +uo{z)m{z) = 
or to Un{z) = b'4)'^_^{z). □ 

By conjugation, we obtain: 

Theorem 2.17. Let Imz > 0. 

(i) Any solution {vn{z)}'^^Q of (12.550 {for n > 2) can be represented as 

v„{z) = p^_^{z)a + qk-i{z)b. (2.59) 

In fact, a = vi{z) and b = —vq[z). 

(ii) A sequence (]2.59p satisfies (j2.55p for all n > 1 if and only ifb = 0. 

(iii) A sequence (j2.59p belongs to l"^ if and only ifvn{z) = ip^_i{z)c for some c G Equivalently, 
a sequence (12.591) belongs to £^ if and only if vi{z) + m{z)vo{z) = 0. 

2.6 Wronskians and the Christoffel— Darboux Formula 

For any two 7W /-valued sequences Vn, define the Wronskian by 

Wn{u, v) = UnAnVn+l " Un+lAl^Vn- (2.60) 

Note that Wn{u,v) = —Wn{v\u'^)'^. If n„(z) and Vn{z) are solutions to (j2.54p and (j2.55p . then 
by a direct calculation, we see that Wn{u{z),v{z)) is independent of n. Put differently, if both 
Un{z) and Vn{z) are solutions to (|2.54p . then Wn{u{z),v{z)'^) is independent of n. Or, if both Un{z) 
and Vn{z) are solutions to (|2.55p . then Wn{u{z)\v{z)) is independent of n. In particular, by a 
direct evaluation for n = 0, we get 

vf^n(p^'-i(^),P^-i(^)^) = iy„(g^_i(z),</^_i(z)t) = 0, 

Wn(pf'-l(^)^pf_l(^)) = I^„(gf'_l(^)^<zf'-l(^)) = o, 

T^n(p^-l(^),9^_l(^)^-) = Wn{p'^_^{z)\q^^^{z)) = 1. 

Let both u{z) and t;(z) be solutions to (|2.54p of the type (|2.58p . namely, 

Un{z) = ap^__^{z) + bq^_^iz), t>„(z) = cp^_^{z) + dq^_^{z). 
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Then the above calculation implies 

Wn{u{z),v{z)^) =aS -bcl 
Theorem 2.18 (CD Formula). For any x,y €z C and n>l, one has 



(x - y) Y.P^{x)p]:{y) = -Wn+M-i{x),p'^_M)- (2-61) 

A;=0 

Proof. Multiplying (j2.30p by Pn{y) on the right and (j2.32p (with y in place of x) by Pn{x) on the 
left and subtracting, we get 

{x - y)p^{x)p^{y) = Wn{p'^^i{x),p':_M) - VF„+i(p!Li(x),pf'_i(y)). 
Summing over n and noting that Wq{p^{x),p^ {y)) = 0, we get the required statement. □ 

2.7 The CD Kernel 

The CD kernel is defined for z, S C by 

n 

Kn{z,w) =Y,Pk{^)Pk{^)^ (2.62) 

n 

= ^pi{-z)^pj^{w). (2.63) 



k=0 



()2.63p follows from (j2.62p and (j2.3ip . Notice that K is independent of the choices o'n,Tn in 
(p:2i|) and that (pleT]) can be written 

{z-w)Kn{z,w) = -Wn+i{p^^iiz),p^_^{w)^). (2.64) 

The independence of Kn of a, r can be understood by noting that if fm is given by ()2.26p . then 

/n 
Kn{z, w) dfi{w)f{w) = ^ p^{z)U (2.65) 

m=0 

SO K is the kernel of the orthogonal projection onto polynomials of degree up to n, and so intrinsic. 
Similarly, if f^^ = {{f,pl;,))L, so 

oo 
m=0 

then, by (12:63]) . 

/n 
f{z)df,{z)Kr,{z,w) = J2 f^f'^PiH- (2-67) 

m=0 

J Kniz,w)df,{w)Kniw,C) = Kniz,C) (2-68) 



One has 
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as can be checked directly and which is an expression of the fact that the map in (j2.65p is a 
projection, and so its own square. 

We will let tt^ be the map of L^{dfj.) to itself given by (or by i^W} ). (I2:6D can then be 

viewed (for z,w € M) as an expression of the integral kernel of the commutator [J, 7r„], which leads 
to another proof of it [175j. 

Let Jn;F be the finite nl x nl matrix obtained from J by taking the top leftmost n2 blocks. It 
is the matrix of TTn-iMxTin-i where is multiplication by x in the {pf}^Zo basis. For y G C and 
7 G C', let ipn,'y{y) be the vector whose components are 

{^nM)3=pf-iiy)7 (2.69) 

for j = 1,2, . . . ,n. We claim that 

[{Jn-F - y)iPna{y)]j = -^jfiAnPnivh (2-70) 

as follows immediately from (I2.32p . 

This is intimately related to (|2.6ip and (|2.64p . For recalling J is the matrix in basis, ^n,'y{y) 
corresponds to the function 

n-l 

J2Pjix)i'Pn,j{y))j~l = Kn{x,y)j. 
3=0 

As we will see in the next two sections, (|2.70p has important consequences. 
2.8 ChristofFel Variational Principle 

There is a matrix version of the Christoffel variational principle (see Nevai |152j for a discussion of 
uses in the scalar case; this matrix case is discussed by Duran-Polo |76j): 

Theorem 2.19. For any non-trivial I x / matrix- valued measure, dfj,, on M, we have that for any 
n, any xq G M, and matrix polynomials Qn{x) of degree at most n with 

Qn{xQ) = 1, (2.71) 

we have that 

{{Qn, Qn))R > Knixo, Xq)"^ (2.72) 

with equality if and only if 

Qn{x) = Kn{x,xo)Kn{xo,xoy^. (2.73) 
Remark. (j2.72p also holds for ((•, but the minimizer is then Kn{xo, xo)~^ Kn{x, xq)- 

Proof. Let Qn^ denote the right-hand side of (|2.73p . Then for any polynomial i?„, of degree at 
most n, we have 

{{Qi'^\Rn))R = Kn{xo,xor^Rn{xo) (2.74) 
because of (I2.65j) . Since Qni^o) = Qn\xo) = 1, we conclude 

{{Qn - Qi?,Qn - Qi?))R = {{Qn,Qn))R " i^„(xo,Xo)-^ (2.75) 

from which ()2.72p is immediate and, given the supposed non- triviality, the uniqueness of minimizer. 

□ 
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With this, one easily gets an extension of a result of Mate-Nevai |146j to MOPRL. (They had 
it for scalar OPUC. For OPUC, it is in Mate-Nevai-Totik [Hlj on [-1, 1] and in Totik [W] for 
general OPRL. The result below can be proven using polynomial mappings a la Totik [188] or Jost 
solutions a la Simon [174 ].) 

Theorem 2.20. Let dfj, be a non-trivial I x I matrix-valued measure on M with compact support, 
E. Let I = (a, b) be an open interval with L C E. Then for Lebesgue a.e. x G L , 

limsup(n + \)Kn{x,x)^'^ < w{x). (2.76) 

Remark. This is intended in terms of expectations in any fixed vector. 

We state this explicitly since we will need it in Section 15.41 but we note that the more detailed 
results of Mate-Nevai-Totik [147], Lubinsky [HI], Simon [T73], and Totik [189] also extend. 

2.9 Zeros 

We next look at zeros of det(P^(z)), which we will prove soon is also det(P^(z)). Following 
[661 I177j . we will identify these zeros with eigenvalues of Jn-F- It is not obvious a priori that 
these zeros are real and, unlike the scalar situation, where the classical arguments on the zeros 
rely on orthogonality, we do not know how to get reality just from that (but see the remark after 
Theorem [225]). 

Lemma 2.21. Let C{z) be an I x I matrix-valued function analytic near z = 0. Let 

k = dim(ker(C(0))). (2.77) 
Then det(C(z)) has a zero at z = of order at least k. 

Remarks. 1. Even in the 1x1 case, where k = 1, the zeros can clearly be of higher order than k 
since cn (z) can have a zero of any order! 

2. The temptation to take products of eigenvalues will lead at best to a complicated proof as 
the cases C{z) = (? g) ^^'^ ^i^) — (i ^q) illustrate. 

Proof. Let ei, . . . , e; be an orthonormal basis with ei, . . . , Cfc G ker(C(0)). By Hadamard's inequal- 
ity (see Bhatia [13]). 

|det(C(z))| < \\C{z)ei\\---\\C{z)ei\\ 
< C\z\'' 

since ||C(2:)ej|| < C|2;| if j = 1, . . . , A; and ||C(z)ej || < d for j = k -\- 1, . . . ,1. □ 

The following goes back at least to [31]; see also [661 [T65l [T771 [T78] . 
Theorem 2.22. We have that 

det0 {P^{z)) = detcni(2 - Jn;F)- (2.78) 
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Proof. By (j2.70p . if 7 is such that Pn{y)l = 0, then ipn^^{y) is an eigenvector for Jn-p with eigenvalue 
y. Conversely, if ip is an eigenvector and 7 is defined as that vector in C' whose components are the 
first / components of then a simple inductive argument shows (p = (pn,-yiy) and then, by ()2.70p 
and the fact that is invertible, we see that Pn{y)l = 0. This shows that for any y, 

dimker(P^(y)) = dimker(J„;F - y). (2.79) 

By Lemma 12.211 if V is any eigenvalue of Jn\F of multiplicity k, then det(P^(z)) has a zero 
of order at least k at y. Now let us consider the polynomials in z on the left and right in (j2.78p . 
Since Jn-^p is Hermitian, the sum of the multiplicities of the zeros on the right is nl. Since the 
polynomial on the left is of degree nl, by a counting argument it has the same set of zeros with 
the same multiplicities as the polynomial on the right. Since both polynomials are monic, they are 
equal. □ 

Corollary 2.23. All the zeros of det{P^{z)) are real. Moreover, 

det{P^{z)) = det{P^{z)). (2.80) 

Proof. Since Jn-p is Hermitian, all its eigenvalues are real, so (|2.78p implies the zeros of det(P^(2;)) 
are real. Thus, since the polynomial is monic. 



det{P^{z)) = det{P^{z)). (2.81) 

By Lemma l2.2f v). we have 

P^{z) = P^{-z)^ (2.82) 
since both sides are analytic and agree if z is real. Thus, 



det{P^{z)) = det(Pi'(z)t) = det{pL{z)) 
proving □ 

The following appeared before in [178 j: see also |165j . 
Corollary 2.24. Let {xn,j}^Li zeros ofdet{Pn{x)) counting multiplicity ordered by 

Xn,l < Xn,2 < ■ ■ < Xn,nl- (2.83) 

Then 

Xn+l,j < Xn,j < Xn+l,j+l. (2.84) 

Remarks. 1. This is interlacing if / = 1 and replaces it for general /. 

2. Using An invertible, one can show the inequalities in (j2.84p are strict. 

Proof. The min-max principle |159j says that 

Xn,j = max min (/, J„;f/)c"!- (2.85) 

dim{L)<j-l 11/11=1 

If P: C("+i)' ^ C"' is the natural projection, then {Pf,Jn+i-FPf)c(n+i)i = {Pf,Jn;FPf)c"i and 
as L runs through all subspaces of C^""*"^^' dimension at most j — 1, P[L] runs through all subspaces 
of dimension at most j — 1 in C"', so (|2.85p implies Xn+ij < Xnj. Using the same argument 
on -Jn;F and -Jn+i-p shows Xj{-Jn-F) > Xj{-Jn+i-F)- But Xj{-Jn-p) = -Xni+i-jiJn;F) and 
Xj{—Jn+i;F) = —X(^n+i)i+i-jiJn+i;F)- That yields the other inequality. □ 
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2.10 Lower Bounds on p and the Stieltjes— Weyl Formula for m 

Next, we want to exploit (|2.7U|) to prove uniform lower bounds on ||p^(y)7|| when y ^ cvh(cj(J)), 
the convex hull of the support of J, and thereby uniform bounds ""^H- We will then use that 

to prove that for z ^ o'(J), we have 



m{z) = lim -p'^{z)-'q'^{z) (2.86) 

n—roo 

the matrix analogue of a formula that spectral theorists associate with Weyl's definition of the 

m- function [193j . although for the discrete case, it goes back at least to Stieltjes |181j . 

We begin by mimicking an argument from |171j . Let H = cv\i{a{J)) = [c — D,c + D] with 

D = \ diam(i7). (2.87) 

By the definition of An, 

Pnll = IK(p1i, (x - c)p^))r\\ < D. (2.88) 

Suppose y ^ H and let 

d = dist{y,H). (2.89) 
By the spectral theorem, for any vector ip £ TLy, 

\{ip,{J-y)^)n^>d\M\ (2.90) 

By (1220]), with if = ipn,^{y), 

\{v,{J-y)v)\ < Pn||||p^(?/)7lll|Pn-i(y)7ll (2.91) 

while 

M' = T.\\pf^y^^\\'- (2.92) 

j=0 

As in [mi Prop. 2.2], we get: 
Theorem 2.25. Ify^H, for any 7, 

\\p'niyh\\>{i){i + {i)T~'^^"\h\\- (2-93) 

In particular, 

\\p^{y)-'\\<i- (2.94) 
Remark. (j2.93p implies det{p^(y)) 7^ if Imy > 0, providing another proof that its zeros are real. 

By the discussion after (j2.52p . if Imz > 0, q^iz) + Pn{z)m{z) is in i"^, so goes to zero. Since 
Pn{z)~^ is bounded, we obtain: 

CoroUary 2.26. For any z G C-|- = {z : Imz > 0}, 

m{z) = lim -p^(z)-ig^(z). (2.95) 

n— >oo 

Remark. This holds for z ^ H. 

Taking adjoints using ()2.82p and m{z)'f = m{z), we see that 

m{z) = lim -q^{z)p^{z)-\ (2.96) 

n~>oo 

p.95p and p.96p are due to [17] , which uses the proof based on the Gauss-Jacobi quadrature 
formula. 
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2.11 Wronskians of Vector- Valued Solutions 

Let a, P be two vector-valued solutions (C') of (|2.55|) for n = 2, 3, . . . . Define their scalar Wronskian 
as (Euclidean inner product on C') 

Wn{a,^) = {an,An(3n+i) " (A„a„+i,/3„) (2.97) 

for 71 = 2, 3, ... . One can obtain two matrix solutions by using a or (3 for one column and for the 
other columns. The scalar Wronskian is just a matrix element of the resulting matrix Wronskian, 
so Wn is constant (as can also be seen by direct calculation). Here is an application: 

Theorem 2.27. Let zq G M\cJess(<^)- For k = 0,1, let he the multiplicity of zq as an eigenvalue 
ofJ^f^K Then, qo + qi<l. 

Proof. If P is an eigenfunction for j(^) and we define /3 by 

Pn = b " = (2.98) 
\Pn-i, n>2, ^ ^ 

then /5 solves ()2.55p for n > 2. If a is an eigenfunction for J = J^^\ it also solves ()2.55p . Since 
On — > 0, /3n ^ 0, and An is bounded, Wn{a,(3) — > as n ^ cxd and so it is identically zero. But 
since /?i = 0, 

= Wi{a,(3) = {ai,AiP2) = {ai,Aj,). (2.99) 
Let V^''^ be the set of values of eigenfunctions of J^^-* at n = 1. (|2.99p says 

C [AiV^'^Y- (2-100) 
Since qk = dim(F(^)) and Ai is non-singular, ()2.100p implies that qo < I — qi- □ 

2.12 The Order of Zeros/Poles of m{z) 

Theorem 2.28. Let zq G M\(Tess(<^)- For A; = 0, 1, let qj. he the multiplicity of zq as an eigenvalue 
of j(^). // qi — qQ> 0, then det(m(2;)) has a zero at z = zq of order qi — qo- If qi — qo < 0, then 
det(m(2;)) has a pole at z = zq of order qo — qi- 

Remarks. 1. To say det(m(z)) has a zero at z = zq of order means it is finite and non-vanishing 
at Zq' 

2. Where dfi is a direct sum of scalar measures, so is m{z), and det(m(2:)) is then a product. 
In the scalar case, m{z) has a pole at zo if J^^^ has zq as eigenvalue and a zero at zo if J^^^ has zq 
as eigenvalue. In the direct sum case, we see there can be cancellations, which helps explain why 
Qi — go occurs. 

3. Formally, one can understand this theorem as follows. Cramer's rule suggests det(r?T,(2;)) = 
det(J^^) — z) / det{J^^^ — z). Even though det(j('^) — z) is not well-defined in the infinite case, we 
expect a cancellation of zeros of order qi and qo- For zq ^ H, the convex hull of cj(j(°)), one can 
use (I2.95P to prove the theorem following this intuition. Spurious zeros in gaps of a{J^^^) make this 
strategy difficult in gaps. 

4. Unlike in Lemma 12.211 we can write m as a product of eigenvalues and analyze that directly 
because m{x) is self-adjoint for x real, which sidesteps some of the problems associated with non- 
trivial Jordan normal forms. 

5. This proof gives another demonstration that qo + Qi < ^- 
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Proof. m{z) has a simple pole at z = zq with a residue which is rank qo and strictly negative 
definite on its range. Let 

f{z) = {z- zo)m{z). 

f is analytic near zq and self-adjoint for z real and near zq. Thus, by eigenvalue perturbation theory 
[1241 1159] . f{z) has / eigenvalues pi{z),. . . , pi{z) analytic near zq with pi, p2, ■ ■ ■ , Pq^ non-zero at zq 
and Pqo+i, ... ,pi zero at zq. 

Thus, m{z) has / eigenvalues near zq, Xj{z) = pj{z)/{z — zq), where Ai,...,Ago have simple 
poles and the others are regular. 

By Proposition 12.141 'm{z)^^ has a simple pole at zq with residue of rank qi (because Ai is 
non-singular), so m{z)~^ has qi eigenvalues with poles. That means qi of Ag^+i, . . . , A; have simple 
zeros at zq and the others are non-zero. Thus, det(m(z)) = 0^=1 ^ji^) has a pole/zero of order 
qo - qi- □ 

2.13 Resolvent of the Jacobi Matrix 

Consider the matrix 



Gnmiz) = {{p^_i, {X - Z) 



Theorem 2.29. One has 



G„^(.) = ^^r^'^^r^^'^' (2.101) 
P„_i(2)V'm-i(^), ifn<m. 



Proof. We have 



^ ^ Gkmiz)Jmn — zGf^n^z), n ^ k, 
m=l 
oo 

y '] JnmGmki^) = ^G'„fc(z), n ^ k. 



m=l 



Fix A; > and let Um{z) = Gkm{z)- Then Um{z) satisfies the equation (j2.54p for n ^ k, and so we 
can use Theorem 12.161 to describe this solution. First suppose k > 1. As Um is an i'^ solution and 
Um satisfies (j2.54p for n = 1, we have 



ak{z)p^_^{z), m<k, 
bkiz)ip^^_-^{z), m>k. 



If A; = 1, (j2.102p also holds true. For m > k, this follows by the same argument, and for m = k = 1, 
this is a trivial statement. Next, similarly, let us consider Vm{z) = Gmk{z)- Then Vm{z) solves 
(j2.55p and so, using Theorem 12.171 we obtain 



IPm_i(^;)cfc(z), mKk, 
ip^_iiz)dkiz), m>k. 



Comparing (|2302|) and (l2T03|) . we find 

akiz)p^_^{z) = Vfc_i(z)ci,n(z), 
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As Pq = Pq = 1, it follows that 

ak{z) = 4^k_^i{z)di{z), Cm{z) = bi{z)i}^,_i{z) 

and so we obtain 



= fTV T/ wr'' ' " (2.104) 



V'n-1 {z)di {z)p^^i (z) if n > m, 
It remains to prove that 

bi{z) = di{z) = 1. (2.105) 
Consider the case m = n = 1. By the definition of the resolvent, 

Gu{z) = {{pliJ-zr'p§U= [^ = miz). 

J X z 

On the other hand, by (|2.104[) . 

Gii(z) = i,^{z)di{z)p^{z) = m{z)di{z), 
Gn{z) = p(^{z)bi{z)i;§{z) =bi{z)m{z), 

which proves (I2.105p . □ 



3 Matrix Orthogonal Polynomials on the Unit Circle 
3.1 Definition of MOPUC 

In this chapter, |U is an Z x ^ matrix-valued measure on dO. ((•, and ((•, are defined as in the 
MOPRL case. Non-triviality is defined as for MOPRL. We will always assume fi is non-trivial. We 
define monic matrix polynomials <I>^, by applying Gram-Schmidt to {1, zl, . . . }, that is, is 
the unique matrix polynomial 2*^1+ lower order with 

{{z''1,^^))r = k = 0,l,...,n-l. (3.1) 

We will define the normalized MOPUC shortly. We will only consider the analogue of what we 
called type 1 for MOPRL because only those appear to be useful. Unlike in the scalar case, the 
monic polynomials do not appear much because it is for the normalized, but not monic, polynomials 
that the left and right Verblunsky coefficients are the same. 



3.2 The Szego Recursion 

Szego |184j included the scalar Szego recursion for the first time. It seems likely that Geronimus 
had it independently shortly after Szego. Not knowing of the connection with this work, Levinson 
|138j rederived the recursion but with matrix coefficients! So the results of this section go back to 
1947. 

For a matrix polynomial P„ of degree n, we define the reversed polynomial P* by 



P:{z) = Z-Pn{l/-zV. 



(3.2) 
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Notice that 

{P*nT = Pn (3.3) 

and for any a S M.i, 

{aPnT = P*a\ {Pno)* = a^P*. (3.4) 

Lemma 3.1. We have 

{{f,9))L = {{9j))l, {{f,9))R = {{9j))^R (3.5) 

and 

{{r,9l)L = {{f,9))i, {{r,9*))R = {{f,9))[. (3.6) 

Proof. The first and second identities follow immediately from the definition. The third identity is 
derived as follows: 

{{f\9*))L= [ e-^5(^)^d^(^)(e"'/(e)1")t 

g{9)Ufi{e)f{e) 

= {{9J))r 
= {{f,9))l 

The proof of the last identity is analogous. □ 



Lemma 3.2. If P-n has degree n and is left- orthogonal with respect to zl, . . . , z"!, then P^ = c(<I>^)* 
for some suitable matrix c. 

Proof. By assumption, 

= {{Pn,zn))L = {{iznr,P:))n = {{z^-n,P:))R for 1 < i < n. 

Thus, P* is right-orthogonal with respect to 1, zl, . . . , and hence it is a right-multiple of 

Consequently, Pn is a left- multiple of (^^)*. □ 

Let us define normalized orthogonal matrix polynomials by 

(/p^ = 99^ = 1, (/.^^ = K^^>^ and ip^ = <^^K^ 

where the k's are defined according to the normalization condition 

along with (a type 1 condition) 

^n+i(^n)"' > and (K^r^K^+i > 0. (3.7) 

Notice that are determined by the normalization condition up to multiplication on the left by 
unitaries; these unitaries can always be uniquely chosen so as to satisfy (|3.7p . 
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Now define 

Pn — '^ni'^n+l) ^ ^^^^ Pn — i'^n+l) ^^ra • 

Notice that as inverses of positives matrices, > and > 0. In particular, we have that 

4 = (Pn-i---Po)"' and = (pI^ . . . p^_^)-\ 
Theorem 3.3 (Szego Recursion), (a) For suitable matrices On^, one has 

^<^n-/'n9^n+l = («n)V^'% (3-8) 

z(^^-<^^+i/>^ = (^^'*(a^)t. (3.9) 

(b) The matrices and are equal and will henceforth be denoted by an- 

(c) = (1 - a]^a„)^/2 pR = (1 _ a„a]^)^/2. 

Proof, (a) The matrix polynomial 299^ has leading term z^+^k^. On the other hand, the matrix 
polynomial p^v^n+i ^^s leading term By definition of p^, these terms are equal. 

Consequently, zy?^ — p!^ip^_^_^ is a matrix polynomial of degree at most n. Notice that it is left- 
orthogonal with respect to zl, . . . , z"l since 

{{zifi: - p^^^+i, z^i))l = {{^t^'''nL - ((p^c^^+i, = 0-0 = 0. 

Now apply Lemma I3.2i The other claim is proved in the same way. 

(b) By part (a) and identities established earlier, 

(a^)t = + (a^)tl 

= {{^^'*,pl:^l:+i))L + {ai:)H{^^,ip^))R 

= ((^^^P^v'^+ik + K)H{^^'\^y))L by m 

= ((^^'*,p^(^^+i + (a^)t^^'*))i 

= {{zipf^, ^Pn'*))R (using the (n + l)-degree *) 

= {{^tlP^ + ^t*{»^)\^^'*fR 

= 0+((^^'*,(^^'*(a^)t))^ 



(c) Using parts (a) and (b), we see that 

1 = {{zipi:,Zip^))L 

{{pnV>n+l + aiv>n'* , PnV>n+l + ai^n'*))l 
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= (Pn)2 + 4((</'n,V'n»iJ«n 

= (/'n)2 + aian- 

A similar calculation yields the other claim. □ 

The matrices an will henceforth be called the Verblunsky coefficients associated with the mea- 
sure d^. Since is invertible, we have 

||q„|| < 1. (3.10) 

We will eventually see (Theorem I3.12p that any set of a^'s obeying (jS.lOp occurs as the set of 
Verblunsky coefficients for a unique non-trivial measure. 

Note that the Szego recursion for the monic orthogonal polynomials is 

z<i>^ - ci>^+i = (K^)-i4(««)t$«'*, 

Z$^-ci>,^^l = cI>^'*(K^)t4(K^)-\ 

so the coefficients in the L and R equations are not equal and depend on all the aj, j = 1, . . . , n. 

Let us write the Szego recursion in matrix form, starting with left-orthogonal polynomials. By 
Theorem 13. 3^ 

'^n+l = [^V'n-V^n'*«]^](/'n)~\ 

= ^(/>n)"Vn " (/>^) "'at (^^'* , (3.12) 

^Si = (P^)-V^'* - z{p^)-'an^l (3.13) 



which implies that 



In other words, 



where 



R,* I — A \<^n, Z) _R,* 



(3.14) 



and = (1 — a^a)"^/^, p^ = {1 — aa^)^^"^. Note that, for z ^ 0, the inverse of A^{a, z) is given by 



A\a,z)-^ = (^ 



:-l(/)-l Z-l(/)-lQt 
ip^Y^a (p^)-i 



which gives rise to the inverse Szego recursion (first emphasized in the scalar and matrix cases by 
Delsarte el al. [37]) 



Matrix Orthogonal Polynomials 



41 



For right-orthogonal polynomials, we find the following matrix formulas. By Theorem 13.31 

^^+1 = z^M-' - (3.15) 
^n+i = ft*{Pnr' - zv^aM)-\ (3.16) 

In other words, 



where 

A^{a,z) 



z{p^)-^ -za{p^y^ 
-at(p«)-i (pi)-i 



For z 7^ 0, the inverse of A^{a, z) is given by 

A^{a,zy^ 
and hence 

= z-V^+i(/^^)-' + z-'^^':,{p^)-'al (3.17) 
3.3 Second Kind Polynomials 

In the scalar case, second kind polynomials go back to Geronimus 1 891 l9T| [92 1 . For n > 1, let us 
introduce the second kind polynomials ipn' by 



V'n(^)= / ^r^J^^{e''>)-^^iz))dpie), (3.19) 



V^^(z) = / dp{e) (^^(e^^) - ip^iz)). (3.20) 

7 6 z 

For n = 0, let us set iPq{z) = iPq{z) = 1. For future reference, let us display the first polynomials 
of each series: 

^liz) = ip(^)-\z - 4), (z) = {z- al){p§)-\ (3.21) 

i^fiz) = {ph)-\z + 4). <(^) = {z + 4)ip^)-\ (3.22) 

We will also need formulas for ' and ipn' , n > 1. These formulas follow directly from the above 
definition and from 



e»^ + l/z\ e^' + z 



e^'^-l/zj e'-'^ - z 



Indeed, we have 



pi6 , 

V'n'*(^) = ^" / rz MO) iv^il/z)^ - vi^ie^y), (3.23) 



ci9 

e*^ + z 



^n'*{z) = / iv^il/z)^ - ^^{e^y) dp{e). (3.24) 



e*^ - z 
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Proposition 3.4. The second kind polynomials obey the recurrence relations 

i^n+ii^) = {pi:)-\zi^i^{z) + at (3.25) 
V'f;*i(^) = {p^rHzani^^^iz) + (3.26) 



and 



V'^+ilz) = (zV«(z) + i;^'*iz)ai){pi^)-\ (3.27) 



i'nU^) = i^t*iz) + zij^{z)anM)-' (3.28) 



for n > 0. 



Proof. 1. Let us check (j3.25p for n > 1. Denote the right-hand side of (j3.25p by il)^_^_i{z). Using 
the recurrence relations for ip^, ipn'* and the definition (j3.19|) of -0^' ^"^^ 



e^^ + z 



Ch_i(^)-C+i(2)= / -^An{e,z)dp{e) 

J e z 

where 

An(^,^)=^^+i(e*')-<^^+i(z) 

- {p^^)-'[z^^^{e^% - z^i{z) + 4z>^(l/z-)t - 4z>^(e^^)t] 

= (p^)-Me^V^(e^') - " ^V'^(^) + 

- z^i{e'') + Z9.^(z) - at ^^'*(z) + a],z>^(e^^)t] 
= - z)^i{e^') + a],(z"e— ^ - l)v^^'*(e^^)]. 

Using the orthogonality relations 

y (^^(e*^)d^(0)e~^™^ = y 99^'*(e^^)(i/i(^)e~*^"'+^^^ = 0, (3.29) 
m = 0, 1, . . . , n — 1, and the formula 



we obtain 



Pn I ^^e^^Ar.{e,z)dm= l{e^'^^{e^')-ay^'*{e^'))dp{e) 



Pn J V>f:+i{e'')dp{e) = 0. 



2. Let us check (|3.26p for n > 1. Denote the right-hand side of (|3.26p by ijj^^i{z). Similarly to 
the argument above, we find 

V^;*i(.)-^f;\(^)= / '^^^B^{e,z)dp{e), 
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where 

Using the orthogonahty relations ()3.29p . we get 

/pid _|_ ^ 
-^Bnie,z)dm = 
6 Z 



3. Relations (j3.25p and (j3.26p can be checked for n = by a direct substitution of (j3.2ip and 
([322]). 

4. We obtain (f3:27p and from (13:251) and (f3:26]l by applying the *-operation. □ 
Writing the above recursion in matrix form, we get 



for left-orthogonal polynomials and 



for right-orthogonal polynomials. 
Equivalently, 



and 

In particular, we see that the second kind polynomials ipn' correspond to Verblunsky coeffi- 
cients {—On}- We have the following Wronskian-type relations: 

Proposition 3.5. For n > and z £ C, we have 

2z^l = (3.31) 

2z"l = + (3.32) 

= ip^{z)iJ^iz)-^P^{z)ip^{z), (3.33) 

= i;^'*iz)^^'*iz) - vp«-*(z)^^'*(z). (3.34) 
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Proof. We prove this by induction. The four identities clearly hold for n = 0. Suppose (j3.3ip - ()3.34p 
hold for some n > 0. Then, 

L I L.* I ; L L,* 

where we used (|3.3ip - (j3.34p for n in the third step. Thus, (j3.3ip holds for n + 1. 
For ([332]) . we note that 

iR,* R , R,* iR 

+ ((p^-* - zan^^){z^lJ^ + ^n'*al)] (Pn)"' 

+ z2a„(V'^(p^ - (p^V'^) - (V^'V^-* - <P^'>^'*)a],] (p^)"' 
= (p^)-i2.'^+i(l - {p^y' = 2z"+il, 

again using p.3ip - (|3.34p for n in the third step. 
Next, 



'Pn+l^n+l - V'n+l'Pn+1 



n "n/ 



= (p^^)-i[(z(p^-4(p^'*)(z^;^ + ^;;'a 

- (zV^ + 4 - v^^'*at )](p«)-i 

= (p^)-i[z2(^^V^ - ^^^Z) - aU^^i^t* - ^n'*^n'*H 

- ZaU<Pn'>n + i^n'^n) + ^(^nV'^'* + 'Pn'*)«]^] (Pn 

= 

which implies first p.33p for n + 1 and then, by applying the *-operation of order 2n + 2, also (j3.34p 
for n + 1. This concludes the proof of the proposition. □ 

3.4 Christoffel— Darboux Formulas 
Proposition 3.6. (a) (CD)-left orthogonal 

n 

(1 - ^z) (p^(e)V^(^) = v^io^^^iz) - ez^^(e)V^(^) 



k=0 
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(b) (CD)-right orthogonal 

n 
k=0 



(c) (Mixed CD)-left orthogonal 

n 

(1 - ^z) ^ ^^{O^ifiiz) = 2 . 1 - - ^z^^iO^^^iz) 



k=0 



2 . 1 - Vn;*i(OVn;*i(^) - ^^+i(6V^+i(^). 



(d) (Mixed CD)-right orthogonal 

(1 - ^z) (^)V'f (0^ = 2 • 1 - <^^'*(z)t/'^'*(o^ - ^"^¥'^(^)V'^(o^ 



fc=o 



Remark. Since the V^'s are themselves MOPUCs, the analogue of (a) and (b), with all t^'s replaced 
by ijj's, holds. 

Proof, (a) Write 



Then, 
and 



Fn+l{^)=A''{0^n,z)F^{z) 



-«(/)-! (p-^)-i y \zip^)-^a -ip^)-^ 

'C^(P^)-' - e^at(p«)-2a -C(p^)-2at + eV(p^)-' 
-za(/)-2 + z{p^)-^a a^p^y^a^ - (p^)-^ 

'C^l 



-1 



J. 



Thus, 



Ft+AO^JF^^^iz) = F^:{i)^A\ar,,O^JA^{ar,,z)F^{z) 



'Hz) 



and hence 



= F„^(OtjX' 
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which shows that the last two expressions in (a) are equal. Denote their common value by Qn{z, 0- 
Then, 

Summing over n completes the proof since Q^i{z,^) = 0. 

(b) The proof is analogous to (a): Write F^{z) = (^(p^{z) (pn'*{z)^. Then, -^^^1(2;) = 
F^{z)A^{an,z) and A^{a, z)JA^{a,0^ = J. Thus, 

F^^,{z)JF^^,0 = F^{z)A^{an,z)JA^{an,O^F^0 
= F^{z)JF^0 

and hence 

^^+1(^)^^+1(0^ - <^n;*l(^)^n;*l(e)^ = ^zcp^{z)^^{0^ - ^t*{z)^n'*{0^ 

which shows that the last two expressions in (b) are equal. Denote their common value by Qn{z, 0- 
Then, 

Q^{z,0 - Qtiiz,0 = (1 - ^z)^^{z)^^0 
and the assertion follows as before. 

(c) Write 

"^-^'^-[-^S'^iz)) 
with the second kind polynomials V'n'^- As in (a), we see that 

F^+iiO^JFn+ii^) = F^iO^A\a^,O^JA\an,z)F^{z) 
= F^iO^JF^iz) 

and hence 

V'^+i(OV^+i(^) + V'^+*i(OV^;V^) = c"^V'^(OV^(^) + v^'*(0V^'*(^). 



Denote 



Then, 



Q^(Z, 0=2-1- Vn4-*l(0Vn;*l(^) - V'n+l(e)Vn+l(^)- 

Qn{z,0 - Qti{z,0 = -V'^'*(OV^(^) - e^V'^(OV^(^) 

+ V'^'*(OV^'*(^) + V'^(OV^(^) 

= (1 - V^(^) 



and the assertion follows as before. 
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(d) Write F^{z) = (ip^iz) - Vn'*(2)) • As in (b), we see that 

F^_,,(z)JK^^,{C)^ = F^(z)A^{a^,z)JA^{a^,O^F^{0^ 
= F^(z)JF^{0^ 

and hence 

With Q^(z,0 = 2 • 1 - (^^;*i(z)V'^;*i(e)t - <^«+i(z)^^+i(e)t, we have 

and we conclude as in (c). □ 



3.5 Zeros of MOPUC 

Our main result in this section is: 

Theorem 3.7. All the zeros of det{ip^{z)) lie in U) = {z : \z\ < 1}. 

We will also prove: 
Theorem 3.8. For each n, 

det{ip^{z)) = det(v9^(z)). (3.35) 

The scalar analogue of Theorem 13.71 has seven proofs in |167j ! The simplest is due to Landau 
|135j and its MOPUC analogue is Theorem 2.13.7 of [167j . There is also a proof in Delsarte et al. 
|37j who attribute the theorem to Whittle |194j . We provide two more proofs here, not only for 
their intrinsic interest: our first proof we need because it depends only on the recursion relation 
(it is related to the proof of Delsarte et al. [37]). The second proof is here since it relates zeros to 
eigenvalues of a cutoff CMV matrix. 

Theorem 3.9. We have 

(i) For z G (9D, all of (fn'*{z), ipn*{z), (p^(^z), fn{z) are invertible. 

(ii) For z G (9B, ip!^{z){ipn'* {zj)~^ and {ipn {z))~^ip^{z) are unitary. 

(iii) For z G D, ipn'*{z) and ipn*{z) are invertible. 

(iv) For z G B, ip^{z){(pn'* (z))^^ and {ipn^ {z))~^(p^{z) are of norm at most 1 and, for n > 1, 
strictly less than 1. 

(v) All zeros of det {(fn'* {z)) and det{ipn* {z)) lie in C \ B. 

(vi) All zeros of det{(p!^{z)) and det{(p!^{z)) lie m B. 

Remark, (vi) is our first proof of Theorem 13. 7i 

Proof. All these results are trivial for n = 0, so we can hope to use an inductive argument. So 
suppose we have the result for n — 1. 

= {ptirHi - ^an-i^ti('^n-*i)-')v=f-*i- (3.36) 
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Since lon-il < if \z\ < 1, each factor on the right of (j3.36p is invertible. This proves (i) and 
(iii) for Lfn'* and a similar argument works for c^n'*- If z = e*^, (p^{e^^) = e*"^93ra'*(e*^)^ is also 
invertible, so we have (i) and (iii) for n. 
Next, we claim that if z G 9B, then 

^^'*(z)t^^.*(z)=^^(z)t^^(z). (3.37) 

This follows from taking z = ^ G (91D) in Proposition I3.6l fa). Given that ipn'*{z) is invertible, this 
implies 

1 = (^^(z)^^'*(z)-i)t(c^^(z)c^«.*(z)-i) (3.38) 

proving the first part of (ii) for n. The second part of (ii) is proven similarly by using Proposi- 
tion [3SIb)._ 
For z e B, let 

F{z) = ^'^{z)^^^*{z)-\ 

Then F is analytic in D, continuous in D, and ||-F(z)|| = 1 on (9B, so (iv) follows from the maximum 
principle. 

Since ipn'*{z) is invertible on B, its det is non-zero there, proving (v). (vi) then follows from 

det((^^(z)) = z"' det((^^'*(l/z)) . (3.39) 

□ 

Let V be the C'-valued functions on 9B and V„ the span of the C'-valued polynomials of degree 
at most n, so 

dim(V„) = C'('^+^). 

Let Voo be the set U„Vn of all C'-valued polynomials. Let 7r„ be the projection onto Vn in the V 
inner product (|1.35p . 
It is easy to see that 

Vn n V^_i = {^^{z)v -.veC'} (3.40) 

since {z\ ^l^{z)v) = for / = 0, . . . , n — 1 and the dimensions on the left and right of ()3.40p coincide. 
Vn n V^-i can also be described as the set of {v'^^^{z))'< for v G C'. 

We define : V„_i — > Vn or Voo Voo as the operator of multiplication by z. 

Theorem 3.10. For all n, we have 

detci(^n(^)) = detv„_i(zl - 7r„_iM,7r„„i). (3.41) 

Remarks. 1. Since H-M^H < 1, (13.411) immediately implies zeros of dQi{(p^{z)) lie in D, and a small 
additional argument proves they lie in D. As we will see, this also implies Theorem 13.81 
2. Of course, 7rn-iM^7r„_i is a cutoff CMV matrix if written in a CMV basis. 

Proof. If Q G V„_fc, then by ([330]), 

Tln-l [{z - zofQ] = {)^{z-zofQ = ^'^{z)v (3.42) 

for some v G Thus writing det(<I>^(2:)) = <^^{z)vi A ■ ■ ■ A ^^{z)vi in a Jordan basis for <I>^(z), 
we see that the order of the zeros of det($^(2;)) at zq is exactly the order of zq as an algebraic 
eigenvalue of iTn-iMziTn-i, that is, the order of zq as a zero of the right side of (|3.4ip . 

Since both sides of ()3.4ip are monic polynomials of degree nl and their zeros including multi- 
plicity are the same, we have proven (j3.4ip . □ 
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Proof of Theorem \3.8[ On the right side of (|3.42p . we can put {^^{z)v^)^ and so conclude (j3.4ip 
holds with $^(2;) on the left. 

This proves ()3.35p if ip is replaced by ^. Since a*aj and aja* are unitarily equivalent, det(/)^) = 
det{p^). Thus, det(K^) = det(K^), and we obtain (|3.35p for ip. □ 

It is a basic fact (Theorem 1.7.5 of |167j ) that for the scalar case, any set of n zeros in D are the 
zeros of a unique OPUC and any monic polynomial with all its zeros in B is a monic OPUC. 
It is easy to see that any set of nl zeros in D is the set of zeros of an OPUC <!>„,, but certainly not 
unique. It is an interesting open question to clarify what matrix monic OPs are monic MOPUCs. 



3.6 Bernstein— Szego Approximation 

Given {aj}^^Q G D", we use Szego recursion to define polynomials (p^-, ipj' for j = 0, 1, . . . , n. We 
define a measure d/iji on dB by 

df^n{e) = [v^{e^')^^{e'y]-'^. (3.43) 
Notice that (|3.37p can be rewritten 

'Pn[e )^n{e y=ipn{e )Vn(e )■ (3.44) 

We use here and below the fact that the proof of Theorem 13.91 only depends on Szego recursion 
and not on the a priori existence of a measure. That theorem also shows the inverse in (I3.43P exists. 
Thus, 

JO 

d/.„(^) = [^^(e^^)ty,^(e*^)]-i-. (3.45) 

Theorem 3.11. The measure d^n is normalized (i.e., /Un(5][D) = 1) and its right MOPUC for 
j = 0,...,n are {y^^j^^Q, and for j > n, 

ipfiz) = z^-^p>l{z). (3.46) 
The Verblunsky coefficients for dfin are 

»j{dl^n) = h' (3.47) 
[0, J >n + l. 

Remarks. 1. In the scalar case, one can multiply by a constant and renormalize, and then prove 
the constant is 1. Because of commutativity issues, we need a different argument here. 

2. Of course, using (fOHjl . ip^ are left MOPUC for d^n- 

3. Our proof owes something to the scalar proof in [80j . 



Proof. Let ((•, be the inner product associated with ^n- By a direct computation, {{ipn, Vn))R 
1, and for j = 0, 1, . . . , n — 1, 



((^^ ^n))R = ^ r e-'^\y^^{e^yr' de 



1 

2^i 



z^~^-\^^'*iz))-Uz = 
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by analyticity of ifn*{z)~^ in B (continuity in D). 

This proves is a MOPUC for d^n (and a similar calculation works for the right side of (j3.46p 
if j > n). By the inverse Szego recursion and induction downwards, {^f}^ZQ are also OPs, and 
by the Szego recursion, they are normalized. In particular, since ip^ = 1 is normalized, d^n is 
normalized. □ 



3.7 Verblunsky's Theorem 

We can now prove the analogue of Favard's theorem for MOPUC; the scalar version is called 
Verblunsky's theorem in [167j after |192j . A history and other proofs can be found in |167j . The 
proof below is essentially the matrix analogue of that of Geronimus ^90j (rediscovered in [371 EO] ) . 
Delsarte et al. |37j presented their proof in the MOPUC case and they seem to have been the first 
with a matrix Verblunsky theorem. One can extend the CMV and the Geronimus theorem proofs 
from the scalar case to get alternate proofs of the theorem below. 

Theorem 3.12 (Verblunsky's Theorem for MOPUC). Any sequence {oij}JLQ E B°° is the sequence 
of Verblunsky coefficients of a unique measure. 

Proof. Uniqueness is easy, since the a's determine the V'j^'s and so the <l>j^'s which determine the 
moments. 

Given a sequence {aj}'jlQ, let dfin be the measures of the last section. By compactness oi I x I 
matrix-valued probability measures on 9B, they have a weak limit. By using limits, {(p^y^Q are 
the right MOPUC for dfj, and they determine the proper Verblunsky coefficients. □ 

3.8 Matrix POPUC 

Analogously to the scalar case (see [151 \W0\ \119\ \172\ I197j ). given any unitary P in M.i, we define 

^^{z-p) = zy.ti{z)-^^':,{z)(3l (3.48) 

As in the scalar case, this is related to the secular determinant of unitary extensions of the cutoff 
CMV matrix. Moreover, 

Theorem 3.13. Fix j3. All the zeros of {ipn{z; f3)) lie on 5B. 
Proof. If \z\ < 1, ip^'*i{z) is invertible and 

VP^(z;/3) = -^^\(z)/3t(l - z/3c^^i*,(z)-Vli(^)) 

is invertible since the last factor differs from 1 by a strict contraction. A similar argument shows 
invertibility if \z\ > 1. Thus, the only zeros of det( • ) lie in 5B. □ 



3.9 Matrix- Valued Caratheodory and Schur Functions 

An analytic matrix-valued function F defined on B is called a (matrix-valued) Caratheodory func- 
tion if -F(O) = 1 and IieF{z) = ^{F{z) + F{z)'^) > for every z G B. The following result can be 
found in Thm. 2.2.2]. 
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Theorem 3.14 (Riesz-Herglotz) . If F is a matrix-valued Caratheodory function, then there exists 
a unique positive semi- definite matrix measure djj, such that 

F{z)= [ MO)- (3.49) 

J e z 

The measure dfi is given by the unique weak limit of the measures dfir{0) = ReF(re*^)^ as r ] 1. 
Moreover, 

oo 
n=l 

where 

cn = I e-'^^'dfiiO). 

Conversely, if dfj, is a positive semi-definite matrix measure, then ()3.49p defines a matrix-valued 
Caratheodory function. 

An analytic matrix- valued function / defined on B is called a (matrix- valued) Schur function if 
f{z)^f{z) < 1 for every z G D. This condition is equivalent to f{z)f{z)^ < 1 for every z € P and 
to 11/(2)11 < 1 for every z G B. By the maximum principle, if / is not constant, the inequalities are 
strict. The following can be found in [167^ Prop. 4.5.3]: 

Proposition 3.15. The association 

fiz) = z-\Fiz)-l)iFiz) + l)-\ (3.50) 
F(z) = (l + z/(z))(l-z/(z))-i (3.51) 

sets up a one-one correspondence between matrix-valued Caratheodory functions and matrix-valued 
Schur functions. 

Proposition 3.16. For z G B, we have 

ReFiz) = (1 - -zfiz)^)-Hl - |z|2/(z)t/(^))(l - ^/(^))-i (3.52) 

and the non-tangential boundary values Re-F(e*^) and /(e*^) exist for Lebesgue almost every 9. 
Write dn{9) = w{6)^ + d^s- Then, for almost every 9, 

w{9) =ReF{e''') (3.53) 

and for a.e. 9, det{w{9)) ^ if and only if fie'^ f{e'^) < 1. 

Proof. The identity (|3.52|) follows from (|3.51|) . The existence of the boundary values of / follows by 
application of the scalar result to the individual entries of /. Then (j3.5ip gives the boundary values 
of F. We also used the following fact: Away from a set of zero Lebesgue measure, det(l — zf{z)) 
has non-zero boundary values by general H°° theory. 

(I3.53P holds for {rj, F{z)r])Qi and {r],diJirf)0 for any 77 G C' by the scalar result. We get (|3.53p 
by polarization. From 

w{9) = (1 - e-'''f{e'y )-\l - fie^yfie^'m - e''f{e''))-' 

it follows immediately that /(e*^) V(e*^) < 1 implies det{w{9)) > 0. Conversely, if /(e*^) V(e*^) < 1 
but not fie^^yfie^^) < 1, then det(l — /(e*^)'l"/(e*^)) = and by our earlier arguments det(l — 
e~^^ f{e^^)^)~^ and det(l — e*^/(e*^))~^ exist and are finite; hence det{vj{9)) = 0. All previous 
statements are true away from suitable sets of zero Lebesgue measure. □ 
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3.10 CoefRcient Stripping, the Schur Algorithm, and Geronimus' Theorem 

The matrix version of Geronimus' theorem goes back at least to the book of Bakonyi-Constantinescu 
[6]. Let F{z) be the matrix- valued Caratheodory function (j3.49p (with the same measure ^ as the 
one used in the definition of ((•, •))_r)- Let us denote 

=V^(z)+V;^(z)F(z), 
uZ{z)=i;^{z)+F{z)^^{z). 

We also define 

n^'*(z)=^^'*(z)-F(z)^^'*(z), 

Proposition 3.17. For any \z\ < 1, the sequences u^{z), u^{z), Un*{z), Un'*{z) are square 
summahle. 



Proof. Denote 

fiO) = , 9{9) 

e ™ — z e-" — ^ 

By the definitions ([339|) ~ (l3:2i|) . we have 

^n(^) = ((/,'/'n»L, 

-uZ^*{z) = z-{{^^,g))R, 
-u^'%z)=z-{{^tg))L, 
u^iz) = {{f,^^))n. 

Using the Bessel inequality and the fact that 1^1 < 1, we obtain the required statements. □ 
Next we will consider sequences defined by 

^"^ - A^{an-i,z) ■ ■ ■ A^iao, z) ( 'f\ (3.54) 



tn J V^O 

where G M.i. Similarly, we will consider the sequences 

{sn, tn) = (so, to)^^(ao, z) ■ ■ ■ A^{an-i,z) (3.55) 

Theorem 3.18. Let z G ID and let f be the Schur function associated with dfj, via (I3.49P and (I3.50p . 
Then: 

(i) A solution of (j3.54p is square summahle if and only if the initial condition is of the form 

= ( ^ 
to) \zf{z)c 

for some matrix c in Mi. 
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(ii) A solution of (j3.55p is square summahle if and only if the initial condition is of the form 

(so,io) = {c,czf{z)) 

for some matrix c. 

Proof. We shall prove (i); the proof of (ii) is similar. 

1. By Proposition 13.171 and (I3.14p . (I3.30p . we have the square summable solution 

(:j = («.t(.) (;°) = . "'im + iy (3.50) 

The matrix d = F{z) + 1 is invertible. Thus, multiplying the above solution on the right by d~^c 
for any given matrix c, we get the "if part of the theorem. 

2. Let us check that (pn'*c is not square summable for any matrix c 7^ 0. By the CD formula 
with ^ = we have 



n-l 



fe=i 

and so 



fe=i 

<^^'*(z)t^^.*(z) > (1 - \z\')^^,{z)^^^,{z) = (1 - \z\^)l. 



Thus, we get 

\\iPn'*c\\R2> (1 - |z|2)Trctc> 0. 
3. Let be any square summable solution to (j3.54p . Let us write this solution as 

'Sn\ [ 'fn<^\ , f i^nb \ Sq + to So - to 

+ > = — - — , b= — - — . (3.57) 



tn J y^Pn aj \-tpn 0/2 2 

Multiplying the solution (j3.56p by b and subtracting from (j3.57p , we get a square summable solution 

ip^iz)ia-Fiz)b) 



(pn*{a - F{z)b) 

It follows that a = F{z)b, which proves the "only if part. □ 

The Schur function / is naturally associated with d/U and hence with the Verblunsky coef- 
ficients ao, ai, a2, ■ ■ ■ ■ The Schur functions obtained by coefficient stripping will be denoted by 
/ij /2, /s, • ■ ■ ) that is, fn corresponds to Verblunsky coefficients On, an+i, «n+2, • • • ■ We also write 
fo^f- 

Theorem 3.19 (Schur Algorithm and Geronimus' Theorem). For the Schur functions /o, /i, f2, ■ ■ ■ 

associated with Verblunsky coefficients ao,ai,a2; ■ ■ ■ , the following relations hold: 

fn+l{z) = Z-\p^r\fn{z) - an] [1 " 4/n(^)]-Vn, (3-58) 

Uz) = {p^)-^[zfn+i{z) + an] [1 + zaifn+i{z)]-^ pt (3-59) 
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Remarks. 1. See (jl.Sip and Theorem 11.41 to understand this result. 
2. p.83p provides an alternate way to write (|3.58|) and (j3.59p . 

Proof. It clearly suffices to consider the case n = 0. Consider the solution of (I3.54p with initial 
condition 

1 



By Theorem 13.181 there exists a matrix c such that 



zfi{z)cj \zfo{z) 

zip^y^ - z{p^y^alfo{z) 
-z{p^r^a, + z{p^)-^h{z))- 

From this we can compute zfi{z): 

zh{z) = [-zip^r'ao + z{p§)-'Mz)] [z{pf;)-' - z{p^)-'alMz)]-' 
= ip^r'ifoiz) - ao] [l-4/o(z)]-Vo 

which is (I338p . 

Similarly, we can express /o in terms of /i. From 



A^{ao,zy^ 



c 



ZMZ)J ^"^'^^ \zfl{zyy 

z-^p^y z-\p(;y'al\ ( c 



z-Hp(i)-'c + {p(^y'alhizy 
{p§y^aoc + {p^yyh{zy 



we find that 



zh{z) = [ip^y'aoc + ip^yyhizy] [z-\p(^y'c + ip(^y'alh{zyy' 
= [ip^r'o.0 + iP^y'zhiz)] [z'\p(^y' + ip{^y'alMz)y' 

= {p^y'[ao + zMz)] [z-H + alMz)y'p(^ 
which gives p.59p . □ 

3.11 The CMV Matrix 

In this section and the next, we discuss CMV matrices for MOPUC. This was discussed first by 
Simon in [170j, which also has the involved history in the scalar case. Most of the results in this 
section appear already in |170j : the results of the next section are new here — they parallel the 
discussion in [1671 Sect. 4.4] where these results first appeared in the scalar case. 
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3.11.1 The CMV basis 

Consider the two sequences Xru G defined by 

X2k{z) = Z-'^if^Ciz), X2k-l{z) = z"^+V2i-l(2)> 

X2k{z) = Z-^Lp^j,{z), X2k-l{z) = Z~^Lp^l^_^{z). 

For an integer /c > 0, let us introduce the fohowing notation: i^ is the (/c + l)th term of the sequence 
0, 1, —1, 2, —2, 3, —3, . . . , and jk is the (A; + l)th term of the sequence 0, —1, 1, —2, 2, —3, 3, . . . . Thus, 
for example, ii = 1, ji = —1. 

We use the right module structure of TL. For a set of functions {/fc(z)}^^Q C its module 
span is the set of all sums ^ fk{z)ak with G Mi. 

Proposition 3.20. (i) For any n > 1, the module span of {xfc}fc=o coincides with the module 
span o/{z*'=}^^Q and the module span of {xk}'^^^ coincides with the module span of {z^^}'^^^. 
(ii) The sequences {Xk}'k=o "'"'^^ {xk}kLo '^'"^ orthonormal: 

{{Xk, Xm))R = {{xk,X.m))R = hm- (3-60) 

Proof, (i) Recall that 

ifniz) = K^z"" + linear combination of{l, . . . , z"^^}, 
(p^'*{z) = (k^)^ + linear combination of{z, . . . , z"}, 

where both and {k^V invertible matrices. It follows that 

Xniz) = 'JnZ^" + linear combination of{z*°, . . . , z*""-^}, 
Xniz) = SnZ-^" + linear combination ofjz-"', . . . , z^"~^}, 

where 7^, Sn are invertible matrices. This proves (i). 
(ii) By the definition of ip^ and 99^, we have 

{{V^,^^))r = {{^t*,^t))R = Sum, (3.61) 
((<p^,z™))jj = 0, m = 0,...,n-l; (((^^•*, = 0, m = 1, . . . , n. (3.62) 

From (j3.61|) with n = m, we get 

{{Xn,Xn))R = {{Xn,Xn))R = 1- 

Considering separately the cases of even and odd n, it is easy to prove that 

((Xn,^™))ij = 0, m = io,ii, . . . ,i„-i, (3.63) 
((x„,z'"))ij = 0, m = jo,ji,...,jn-i- (3.64) 

For example, for n = 2k, m £ {zq, . . . , i2k-i} we have m + k £ {1,2, . . . , 2k} and so, by (j3.62p . 

{{Xn, Z^))r = iz-WC^ Z^))r = i^\C^Z^^^))K = 0. 

The other three cases are considered similarly. From ()3.63p . (I3.64|) . and (i), we get ()3.60p for 
k ^ m. □ 
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From the above proposition and the ||-||oo-density of Laurent polynomials in C((9B), it follows 
that {xfclfe^o {^fclfc^o right orthonormal modula bases in TL, that is, any element f £ Ti. 
can be represented as a 

oo oo 

/ = X] XkiiXk, f))R = ^ Xk{{xk, f))R. (3.65) 

k=0 k=0 

3.11.2 The CMV matrix 

Consider the matrix of the right homomorphism f{z) ^ zf{z) with respect to the basis {Xk}- 
Denote Cnm = {{Xn, zXm))R- The matrix C is unitary in the following sense: 

oo oo 

^kn^km = ^^^nkCmk ~ ^nm^- 

k=0 k=0 

The proof follows from (|3.65p : 



^nml = {{zXn,ZXm))R = (T^XkiiXk, ZXn)) R, ZXm^ = "^Cln^. 

^k=0 R k=0 

6nm^ = {{zXn, ZXm)) R = \ ^ XkiiXk, ZXn)) R, ZXm ^ = ^ C„fcC, 



We note an immediate consequence: 
Lemma 3.21. Let \z\ < 1. Then, for every m > 0, 

oo oo 

'^Xn{z)Cnm = ZXm{z), ^ Cm„Xn (1/^)"^ = ZXm{'^/z)'^- 

n=0 n=0 

Proof. First note that the above series contains only finitely many non-zero terms. Expanding 
f{z) = ZXn according to (I3.65p . we see that 

oo oo 

ZXn{z) = J2^k{z){{Xk,ZXn))R = ^Xfel^^jC^ 
fc=0 k=0 

which is the first identity. Next, taking adjoints, we get 



-ZXniz)^ = J24nXk{z)^ 



k=0 



which yields 



oo 



Z^CmnXn{z)^ = ^Cmn^Cl^Xk{z)^ 



n=0 n=0 k=0 

oo / oo 



oo y (JXJ \ 

E(EC-nCl)xk{z)' 

oo 

^SmkXkiz)'' = Xm{z)^- 



k=0 

Replacing z by we get the required statement. □ 
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3.11.3 The £A^-representation 

Using (|3.65|) for / = Xmi we obtain: 

oo oo 
Cnm = {{Xn,ZXm))R = '^{{Xn, ZXk)) R{{xk, Xm)) R = ^nkMkm- 

Denote by G(a) the 21 x 21 unitary matrix 

0(a) = \ ^R 



(3.66) 



k=0 



-a 



Using the Szego recursion formulas (j3.15p and (I3.18p . we get 

Z^n = fn+lPn + fn*(^L 

L,* L.* L R 

"Pn+l ='Pn Pn- V^n+l^n- 



X2fca2fc + X2k+lP2k^ 



Taking n = 2k and multiplying hy z ^, we get 

ZX2k 

ZX2k+l = X2kP2k - X2k+ia2k- 

It follows that the matrix C has the structure 

c = e(ao) e 9(02) e 6(04) e • • • . 

Taking n = 2A; — 1 in (j3.67p . (j3.68p and multiplying by z~^, we get 

X2k-1 = X2k~ialf,^i + X2kP2k-l^ 
X2k = X2k-lP2k-l ~ X2kOi2k-l- 

It follows that the matrix has the structure 

M = le G(ai) e e(a3) © • • • . 
Substituting this into ()3.66p . we obtain: 



Po 








— aoQ;| 
t R 
(^2Pl 

P2P1 





PoPi 








R t 
alp3 










P2P3 

-02/93 

-a\a3 



V 



(3.67) 
(3.68) 



(3.69) 



(3.70) 



We note that the analogous formula to this in |170j . namely, (4.30), is incorrect! The order of the 
factors below the diagonal is wrong there. 
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3.12 The Resolvent of the CMV Matrix 

We begin by studying solutions to the equations 



oo 

CmkWk 

k=0 

oo 

m 

k=0 



ZWr. 



ZWr, 



m > 2, 



m > 1. 



(3.71) 
(3.72) 



Let us introduce the following functions: 



Xn[z) 

T2n-l(^;) 

y2n(z) 
y2n-l{z) 
Pn{z) 

vr„(z) 



Xni^/z) 



-T,„.(l/z)t 



z ->^„„i(z), 



^2n{i^/Z)' = Z 

-T2„-l(l/z)^ = -z 
yn{z) + Xn{z)F{z), 

Tn{z) + F{z)xn{z). 



V2l;-i(^), 



Proposition 3.22. Lei z G B \ {0}. 

(i) For each n > 0, a pair of values {w2mW2n+i) uniquely determines a solution Wn 
Also, for any pair of values {w2n,u!2n+i) in -Mi, there exists a solution Wn to (I3.72p 
values at {2n, 2ra + 1). 

(ii) The set of solutions Wn to (13.720 coincides with the set of sequences 

Wn{z) = aXn{z) + b-Kn{z) 

where a,b range over Mi. 

(iii) A solution (I3.73|) is in l"^ if and only if a = 0. 

(iv) A solution ()3.73p obeys ()3.72p for all m> if and only if b = 0. 

Proposition 3.23. Let z £lD)\ {0}. 

(i) For each n > 1, a pair of values {w2n~i,W2n) uniquely determines a solution Wn 
Also, for any pair of values {w2n~i,'W2n) in -Mi, there exists a solution Wn to ()3.7ip 
values at (2n — 1, 2n). 

(ii) The set of solutions Wn to (|3.7ip coincides with the set of sequences 

Wniz) = Xniz)a + Pniz)b 

where a, b range over A4i. 

(iii) A solution ()3.74p is in i'^ if and only if a = 0. 

(iv) A solution (I3.74p obeys (j3.7ip for all m > if and only if b = 0. 

Proof of Proposition \3.2^ (i) The matrix C — z can be written in the form 

fAo Bo •••\ 
Ai Bi ■ 
^2 ^2 • 



to IK72h . 

with these 



(3.73) 



to (HjH). 

with these 



(3.74) 



V ■■ 



(3.75) 
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where 



Bn 



P2n'^2n+1 P2nP2n+l 



Define Wri = {w2n-,W2n+\) for n = 0, 1, 2, ... . Then (|3.72|) for m = 2n + 1, 2n + 2 is equivalent to 

It remains to prove that the 21 x 21 matrices Aj, Bj are invertible. Suppose that for some x, y E C^, 
An (y) = (q) . This is equivalent to the system 

<^2nP2n-lX - "L^an-iy - Zy = 0, 
P2nP2n-l^ - P^n«2n-iy = 0. 

The second equation of this system yields P2n-i^ — '^2n~iy (since is invertible), and upon 
substitution into the first equation, we get y = x = 0. Thus, kei{An) = {0}. In a similar way, one 
proves that keic{B„) = {0}. 

(ii) First note that Wn = Xn is a solution to (j3.72p by Lemma [3.211 Let us check that Wn = T„ 
is also a solution. If Ukm = {—^)^^km, then {UCU)km for m > 1 coincides with the CMV matrix 
corresponding to the coefficients {— q„}. Recall that ipn^ are the orthogonal polynomials (fn'^, 
corresponding to the coefficients {— a„}. Taking into account the minus signs in the definition of 
T„, we see that Wn = T„ solves (j3.72p for m > 1. It follows that any Wn of the form (j3.73p is a 
solution to (13721) . 

Let us check that any solution to (|3.72p can be represented as ()3.73p . By (i), it suffices to show 
that for any wq, wi, there exist a,b £ Aii such that 

aXo{z) + fevro(z) = wq, 
axiiz) + bTTi{z) = wi. 

Recalling that xo = li '^o = —1) '^li^) = {z + aQ)(/9^)^^, Xii^) = — Oio){Po )~^ ^ see that the 
above system can be easily solved for a,b if z ^ 0. 

(iii) Let us prove that the solution 7r„ is square integrable. We will consider separately the 
sequences 7r2n and 7r2n-i and prove that they both belong to By (|3.20p and (|3.23p . we have 



^P^iz) + F{z)^^{z) = I 4^ dm^Zie''), (3-76) 



z 

- nz)^n'*i^) = / df,i9)^i:{e^y. (3.77) 



Taking n = 2/c in (I3.77P and n = 2fc — 1 in (I3.76p . we get 

7r2kiz) = z^ I 4^ dm^^2k{e'')\ (3.78) 
J a z 
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/pid _|_ 2 
-^d^l{e)^!^,_,[e''). (3.79) 
e z 

As ip2k is an orthonormal sequence, using the Bessel inequality, from (|3.78p we immediately get 
that ■K2k is in 

Consider the odd terms TT2k-i- We claim that 

z-'^' [ 4^ df,{9)^^,_,{e^') = [ ^ d^(0)e^(-'=+i)Vf,-i(e^^). (3.80) 
J e z J a z 

Indeed, using the right orthogonality of 'f2k-i e*"^^, m = 0, 1, . . . , 2A; — 2, we get 



I oo 

2 z'^e'^'^^l^i 



m imB R 
m=l II R 



and 



/ 



m=2k-l 



R 



m=l ^ II R 



2 ^ z-e^-^(^f,_l 

m=2A:-l ^ 

which proves dSSQ]). The identities (jlH)]) and (ISlTQl) yield 

and, since X2k-i is a right orthogonal sequence, the Bessel inequality ensures that ■K2k-i{z) is in 
f. Thus, 7rfc(z) is in i"^. 

Next, as in the proof of Theorem 13. 18^ using the CD formula, we check that the sequence 
Ilv5n'*(-z)||ij is bounded below and therefore the sequence X2n{z) is not in This proves the 
statement (iii). 

(iv) By Lemma |3.21|, the solution Xn{z) obeys ()3.72p for all m > 0. It is easy to check directly 
that the solution iTniz) does not obey (|3.72p for m = if z 7^ 0. This proves the required 
statement. □ 

Proof of Proposition \3 .23[. (i) For j = 1,2,..., define Wj = {w2j-i-,W2j)- Then, using the block 
structure (j3.75p . we can rewrite (j3.7ip for m = 2j, 2j + 1 as AjWj + BjWj+i = 0. By the proof of 
Proposition 13.22] the matrices Aj and Bj are invertible, which proves (i). 

(ii) Lemma 13.211 ensures that Xn{z) is a solution of (j3.7ip . As in the proof of Proposition 13.221 
by considering the matrix {UCU)km, one checks that yn{z) is also a solution to (j3.7ip . 
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Let us prove that any solution to (|3.7ip can be represented in the form (j3.74p . By (i), it suffices 
to show that for any wi, W2, there exist a,b £ Mi such that 

xi{z)a +pi{z)b = wi, 

X2{z)a +P2{z)b = W2. 

We claim that this system of equations can be solved for a, b. Here are the main steps. Substituting 
the definitions of Xniz) and Pn{z), we rewrite this system as 

(Pi '*{a + F{z)b) - '4)^'*{z)b = zwi, 
ip^{a + F{z)b) + '>p2{z)b = zw2. 

Using Szego recurrence, we can substitute the expressions for ip2, ip2^ which helps rewrite our 
system as 

(/3f'*(a + F{z)b) - '4}i'*{z)b = zwi, 

iff {a + F{z)b) + Tpi{z)b = p[w2 + a\wi. 

Substituting explicit formulas for 99^'*, (/^f, V'l') and expressing F{z) in terms of f{z), we can 

rewrite this as 

{p^rHl - aoz)a + 2z{p^)-\f{z) - qo)(1 - zf{z)y^b = zwi, 
{p^r\z - a|))a + 2z{p^r\l - 4/(z))(l - zf{z))~^b = p[w2 + a\wi. 

Denote 

oi = (/Oo^)~^(l - aQz)a, 

61 = 2z{p^r\l - alf{z)){l - zf{z))-^b. 

Then in terms of ai, 61, our system can be rewritten as 

ai + Xibi = zwi, 

X2ai + 61 = P1W2 + a{wi, 

where 

Xi = {p^)-\f{z) - ao)(l - ao/(^))-V^, 
X2 = {pl;)-\z-al){l-aozr^p^. 

Since ||/(^;)|| < 1 and \z\ < 1, we can apply Corollary [L5l which yields ||Xi|| < 1 and IIX2II < 1. It 
follows that our system can be solved for ai, 61. 

(iii) As Pn{z) = — 7r„(l/z)^, by Proposition 13.221 we get that pn{z) is in l"^. In the same way, as 
Xn{z) = Xni^/z)^ , we get that is not in l"^. 

(iv) By Lemma I3.2H the solution Xn{z) obeys ()3.7ip for all m > 0. Using the explicit formula 
for Uniz), one easily checks that the solution yn{z) does not obey (I3.7ip for m = 0, 1. □ 
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Theorem 3.24. We have for z e B, 



J (2z) Xk{z)'Ki{z), I > k or k = I even, 



{2z) pk{z)xi{z), k > I or k = I odd. 



Proof. Fix z e D. Write Gk,i{z) = [{C — z)~^]k,i. Then G.j{z) is equal to (C — z)~^6i, which means 
that Gk,i{z) solves (j3.7ip for m ^ I. Since G.^i{z) is l"^ at infinity and obeys the equation at m = 0, 
we see that it is a right-multiple of p for large k and a right-multiple of x for small k. Thus, 



Gk,i{^) 

Similarly, 

Equating the two expressions, we find 



Xk{z)ai{z), k < I or k = I even, 

Pk{z)bi{z), k > I or k = I odd. 

bk{z)iTi{z), k < I ox k = I even, 

dk{z)xi{z), k > I or k = I odd. 



Xk{z)ai{z) = bk{z)Tri{z) k < I or k = I even, (3.81) 

Pk{z)bi{z) = ai:{z)xi{z) k > I or k = I odd. (3.82) 

Putting = in (j3.8ip and setting bo{z) = ci{z), we find ai{z) = ci{z)-ki{z). Putting / = in 
(|3.82p and setting bo{z) = C2{z), we find Pk{z)c2{z) = dk{z). Thus, 



Gk,ii^) 



Xk{z)ci{z)7ri{z), k < I or k = I even, 
Pk{z)c2{z)xi{z), k > I or k = I odd. 



We claim that ci{z) = C2{z) = {2z) ^1. Consider the case k = I = 0. Then, on the one hand, by 
the definition, 



e^^ + z 



e*^ - z 



dn{e 



10 \ 



(2z)-i 

= {2z)-\F{z) - 1) (3.83) 

and on the other hand, 

Gofiiz) = xo{z)ci{z)tto{z) = ci{z){F{z) - 1). 

This shows ci(z) = {2z)~^l. Next, consider the case A; = 1, / = 0. Then, on the one hand, by the 
definition, 

Gi,o{z) = {{xi,{e''-zr\o))R 

and on the other hand, 

Gi,o{z) =pi{z)c2{z)xo{z). 
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A- 


= ip§)-' ^ 


A- 


= iP^r' ' 


1 



Let us calculate the expressions on the right-hand side. We have 

Piiz)c2iz)xo{z) = ip§)-\-z-' -ao + (z-^ - ao)Fiz))c2iz) (3.84) 

and 

{{xUe''-z)-'xo))R = 

' - z)-i - z-^e-'^ - ao(e^^ - z)-^] df,{9) 

{F{z) - 1) - ^ MF{z) - 1) - ^ / dm 
Taking into account the identity 

j e-'^^dnie) = ao 
(which can be obtained, e.g., by expanding {{ipi, y^Q))R = 0), we get 

{{Xi, ie'' - z)-\o))r = ^ {p^)-\-z-^ - ao + {z^^ - a^)F{z)). 

Comparing this with (|3.84p . we get 02(2) = {2z)^^l. □ 

As an immediate corollary, evaluating the kernel on the diagonal for even and odd indices, we 
obtain the formulas 

^2n{e'') y^Lie^y = vL{z)nt*{z), (3.85) 

^^l-l{e'y^^^l-Ae'') = -:^ut.^{z)^^^^^ (3.86) 

Combining this with (j3.3ip and (j3.32p . we find 

ui{z)^l^\z) + ^i{z)ui^\z) = 2z", (3.87) 
^l^\z)u^{z) + ul^\z)^fX^) = 2z\ (3.88) 

3.13 Khrushchev Theory 

Among the deepest and most elegant methods in OPUC are those of Khrushchev [1251 1126} 1101] . 
We have not been able to extend them to MOPUC! We regard their extension as an important 
open question; we present the first very partial steps here. 
Let 

O = {0 : det?i;(e) > 0}. 

Theorem 3.25. For every n > 0, 

{Q : fnie^yfnie'') < 1} = n 

up to a set of zero Lebesgue measure. 
Consequently, 

.. „ , .. de \n\ 
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Proof. Recall that, by Proposition 13.161 up to a set of zero Lebesgue measure, 

{0 : /o(e^')Vo(e^') < 1} = : detwiO) > 0} 

so, by induction, it suffices to show that, up to a set of zero Lebesgue measure, 

{0 ■■ foie'yUe'') < 1} = : hie'^ Me'') < 1}. 

This in turn follows from the fact that the Schur algorithm, which relates the two functions, 
preserves the property g^g < 1. 

Notice that away from i7, fn{e^^) has norm one and therefore, 

i8,.. de ^ f „ „ , d9 f ^ dO 



which yields (fM]) . □ 
Define 

bn{z;dfi) = ip^{z;dfi)(p^'*{z;dfi)~'^. 

Proposition 3.26. (a) bn+i = - ai){l - zanhn)~^Pn- 

(b) The Verhlunsky coefficients ofhn are {—a\_i, —a\_2: • • • > ~Oq, 1). 

Proof, (a) By the Szego recursion, we have that 

1 Li R,* \ — 1 

= {p^)~\zip^ - ai^y){ipy - zan^^-^pl 

= ipi^r'iz^^Av^'*)-' - ay^'*{^^n'') 

= {Pnr\zbn - al){l - Zanbn)~^Pn- 

(b) It follows from part (a) that the first Verblunsky coefficient of 6„ is —Ci\^_i and that its first 
Schur iterate is ftn-ii compare Theorem 13.191 This gives the claim by induction and the fact that 
6o = 1. □ 



4 The Szego Mapping and the Geronimus Relations 

In this chapter, we present the matrix analogue of the Szego mapping and the resulting Geron- 
imus relations. This establishes a correspondence between certain matrix-valued measures on the 
unit circle and matrix-valued measures on the interval [—2, 2] and, consequently, a correspondence 
between Verblunsky coefficients and Jacobi parameters. Throughout this chapter, we will denote 
measures on the circle by dpic and measures on the interval by dpj. 

The scalar versions of these objects are due to Szego |182j and Geronimus [90j. There are four 
proofs that we know of: the original argument of Geronimus [90] based on Szego's formula in |182j . 
a proof of Damanik-Killip [32] using Schur functions, a proof of Killip-Nenciu |127j using CMV 
matrices, and a proof of Faybusovich-Gekhtman [81j using canonical moments. 
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The matrix version of these objects was studied by Yakhlef-Marcellan [199| who proved The- 
orem 14.21 below using the Geronimus-Szego approach. Our proof uses the Kilhp-Nenciu-CMV 
approach. In comparing our formula with |199j . one needs the following dictionary (their objects 
on the left of the equal sign and ours on the right): 

Hn = -"i+i, 

En = Bn+1- 

Dette-Studden [33] have extended the theory of canonical moments from OPRL to MOPRL. 
It would be illuminating to use this to extend the proof that Faybusovich-Gekhtman gave of 
Geronimus relations for scalar OPUC to MOPUC. 

Suppose dfiQ is a non-trivial positive semi-definite Hermitian matrix measure on the unit circle 
that is invariant under 6 ^ —9 (i.e., z z = z~^). Then we define the measure dfij on the interval 
[-2,2] by 

I f{x)dfxj{x) = I f{2cose)dfic{e) 

for / measurable on [—2,2]. The map 

Sz : dfic ^ dfij 

is called the Szego mapping. 

The Szego mapping can be inverted as follows. Suppose dfij is a non-degenerate positive semi- 
definite matrix measure on [—2,2]. Then we define the measure dfic on the unit circle which is 
invariant under 9 —9 hy 

J g{9)dfic{0) = j 5 (arccos(x/2)) (i/i/(x) 

for g measurable on with g{9) = g{—9). 

We first show that for the measures on the circle of interest in this section, the Verblunsky 
coefficients are always Hermitian. 

Lemma 4.1. Suppose d^c "is a non-trivial positive semi-definite Hermitian matrix measure on the 
unit circle. Denote the associated Verblunsky coefficients by {a„}. Then, dfic "is invariant under 
9 —9 if and only if a\i = c^n for every n. 

Proof. For a polynomial P, denote P{z) = Piz)^ ■ 

1. Suppose that dfic is invariant under 9 —9. Then we have 

{{f,9))L = {{gJ~))R 

for all /, g. Inspecting the orthogonality conditions which define and we see that 

1.^ = $^ and {{<^t ^nk = P^, ^^))r- (4.1) 



Next, we claim that 



(4.2) 
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Indeed, recall the definition of k„ : 

i^n = Un{{^t ^n))L^^^ is Unitary, > 0, 4 = 1> 

= p^, is unitary, {k^)~^k^+i > 0, = 1. 

Using this definition and (j4.ip . one can easily prove by induction that Vn = Un and therefore ()4.2p 
holds true. 

Next, taking z = in (|3.1ip . we get 

a. = -(K«)-icI.^^+i(0)t(K^)t, 

From here and ()4.ip . ()4.2p . we get a„ = On. 

2. Assume On = for all n. Then, by Theorem I3.3r c). we have = p^. It follows that in 
this case the Szego recurrence relation is invariant with respect to the change t/?^ i— > (p^, ip^ i— > (pl{. 
It follows that ifl^ = if^, ip^ = P>n- In particular, we get 

{{Vn,v';n))L = {{^t^n))R- (4-3) 

Now let / and g be any polynomials; we have 

n n 

and a similar expansion for g. Using these expansions and ()4.3p . we get 

{{f,9))L = {{g, f))R for all polynomials /, g. 

From here it follows that the measure dfj-c is invariant under 6 —9. □ 

Now consider two measures dfic and dfij = Sz{dp,c) a-nd the associated CMV and Jacobi 
matrices. What are the relations between the parameters of these matrices? 

Theorem 4.2. Given dfic o-i^d djjLi = Sz{d^c) o-s above, the coefficients of the associated CMV 
and Jacobi matrices satisfy the Geronimus relations: 

Bk+l = \/l - a2k-l Ct2k \/l - Oi2k-l - \/l + Ol2k-l 0!2k-2 \/l + a2fc-l , (4.4) 



Ak+1 = \/l - a2k-i y 1 - a^fc \/l + "2fe+i ■ (4.5) 

Remarks. 1. For these formulas to hold for A; = 0, we set a_i = —1. 

2. There are several proofs of the Geronimus relations in the scalar case. We follow the proof 
given by Killip and Nenciu in [127] . 

3. These A^s are, in general, not type 1 or 2 or 3. 
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Proof. For a Hermitian I x I matrix a with ||a|| < 1, define the unitary 21 x 21 matrix S{a) by 

Of \ _ ^ f - + a 

Since = a, the associated and coincide and will be denoted by p. We therefore have 

e(a) = f« P 

\p -a 

and hence, by a straightforward calculation, 



^ ^ ^ ^ ^ 2 VVl + a VI - a y VP -ay V-Vl + « VI - a 

-1 
1 

5 = 1 e 5(ai) e Sias) 

7^ = 10 (-1)0 10 (-1)6 



Thus, if we define 

and 

it follows that (see IKm\\ ) 



SMS^ = n. 

The matrix C^A + is unitarily equivalent to 

A = S{CM + MC)S^ = SCS^TZ + TZSCS\ 

Observe that A is the direct sum of two block Jacobi matrices. Indeed, it follows quickly from 
the explicit form of TZ that the even-odd and odd-even block entries of A vanish. Consequently, A 
is the direct sum of its odd-odd and even-even block entries. We will call these two block Jacobi 
matrices J and J, respectively. 

Consider C as an operator on Ti^. Then dpc is the spectral measure of C in the following sense: 

[C"^]o,o = j e'^^'dpciO); 

see ()3.83p . Then, by the spectral theorem, (i/ic(— ^) is the spectral measure of and so dpi is 
the spectral measme of C +C-'^ = CM + {CM)-^ = CM +MC. Since S leaves {1 •••)* 
invariant, we see that dpi is the spectral measure of J. 

To determine the block entries of J, we only need to compute the odd-odd block entries of A. 
For > 0, we have 

A I n—, n A (-ot2k-2 ^ \ ( VI + a2fc-A 

A2k+l,2k+l = (VI + «2A:-1 VI - "2fc-lj n /I 

= Vl - «2fc-l a2fc\/l ~ Oi2k-l - Vl + a2A:-l a2fc-2\/l + «2fe-l 
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and 

0\ f^l + a2k+? 

P2k oy \vi - oi2k+i) 

= Vl - Oi2k-l\J'^- O^lkV^ + «2fc+l ■ 

The result follows. □ 
As in \127\ 1168] . one can also use this to get Geronimus relations for the second Szego map. 

5 Regular MOPRL 

The theory of regular (scalar) OPs was developed by Stahl-Totik |180j generalizing a definition 
of Ullman [190j for [—1,1]. (See Simon [171j for a review and more about the history.) Here we 
develop the basics for MOPRL; it is not hard to do the same for MOPUC. 

5.1 Upper Bound and Definition 

Theorem 5.1. Let djjL he a non-trivial I x / matrix-valued measure on M with E = supp((i/^) 
compact. Then [with C{E) = logarithmic capacity of E) 

limsup |det(^i • • • ^n)|^/" < C[E)K (5.1) 

n—>oo 

Remarks. 1. |det(^i • • ■ An)\ is constant over equivalent Jacobi parameters. 

2. For the scalar case, this is a result of Widom |195j (it might be older) whose proof extends 
to the matrix case. 



Proof. Let T„ be the Chebyshev polynomials for E (see [171^ Appendix B] for a definition) and let 
Tn'^ be r„ (8) 1, that is, the I x 
monic so, by (f2J2]l and ^2^, 



Tn'^ be r„ (8) 1, that is, the I x I matrix polynomial obtained by multiplying Tn{x) by 1. Tn^ is 



l/2n 

|det(Ai...A„)|i/"< 



det 



< sup \ Tn{x' 



il/n 



By a theorem of Szego [l83], sup„|r„(x)|i/" ^ C{E), so dST]) follows. □ 

Definition. Let dfi be a non-trivial / x / matrix-valued measure with E = supp{d^i) compact. We 
say is regular if equality holds in (jS.ip . 



5.2 Density of Zeros 

The following is a simple extension of the scalar results (see [71] or |17H Sect. 2]): 
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Theorem 5.2. Let dfi be a regular measure with E = supp((i/i). Let dvn he the zero counting 
measure for det{p^{x)), that is, if {xj^^}^Li are the zeros of this determinant {counting degenerate 
zeros multiply), then 

nl 



dv, 



Then dvn converges weakly to dpE, the equilibrium measure for E. 

Remark. For a discussion of capacity, equilibrium measure, quasi-every (q.e.), etc., see [171^ Ap- 
pendix A] and [Ha [TMfTSH] . 



Proof. By (fZSO]) and ([IJS]), 



|det(p^(x))|>(:^)'(l + (^)Y"-^)'/^ (5.3) 



so, in particular, 

lim^inf |det(p^(x))|i/"' > (l + {iff^ > 1. (5.4) 

But 

|det(p^(x))| = |det(^i • • • ^„)r' exp(-n/$,„(x)) (5.5) 



where is the potential of the measure v. Let z/qo be a limit point of f„ and use equality in (j5.ip 
and (15. 4p to conclude, for x ^ cvh(S), 

exp(-f.^(x))>C(i?) 

which, as in the proof of Theorem 2.4 of [171], implies that Voo = Pe- D 

The analogue of the almost converse of this last theorem has an extra subtlety relative to the 
scalar case: 

Theorem 5.3. Let dfj, be a non-trivial I x / matrix-valued measure on M with E = supp{dfi) 
compact. If dvn dpE, then either p is regular or else, with dp = M{x) dpt^ix), there is a set S 
of capacity zero, so det(M(2;)) = for dptr-a.e. x ^ S. 

Remark. By taking direct sums of regular point measures, it is easy to find regular measures where 
det(M(x)) = for d/itr-a.e. x. 

Proof. For a.e. x with det(M(x)) / 0, we have (see Lemma [52] below) 

P^(x) <C(n + l)l. 

The theorem then follows from the proof of Theorem 2.5 of |171] . □ 
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5.3 General Asymptotics 

The following generalizes Theorem 1.10 of |171| from OPRL to MOPRL — it is the matrix analogue 
of a basic result of Stahl-Totik [180| . Its proof is essentially the same as in jl71j . By cJess(M)) '^^ 
mean the essential spectrum of the block Jacobi matrix associated to /i. 

Theorem 5.4. Let E C M be compact and let fi be an I x I matrix-valued measure of compact 
support with (Tess(/^) = E. Then the following are equivalent: 

(i) fi is regular, that is, limn,^oo|det(Ai • • ■ An)\^^"' = C{E)K 

(ii) For all z in C, uniformly on compacts, 

limsup|det(p^(z))|^/" < e^^^^\ (5.6) 

(iii) For q.e. z in E, we have 

limsup|det(p^(z))|i/" < 1. (5.7) 

Moreover, if (i)-(iii) hold, then 

(iv) For every z G C \ cvh(supp(d/i)), we have 

lim |det(p^(z))|i/" = e^^(^). (5.8) 

n^oo 

(v) For q.e. z € dO-, 

lim sup |det(p^(z))|^/" = 1. (5.9) 

n— »oo 

Remarks. 1. Ge, the potential theorists' Green's function for E, is defined by Ge{z) = 
-\og{C{E))-^p^{z). 

2. There is missing here one condition from Theorem 1.10 of |171j involving general polynomials. 
Since the determinant of a sum can be much larger than the sum of the determinants, it is not 
obvious how to extend this result. 

5.4 Weak Convergence of the CD Kernel and Consequences 

The results of Simon in [174] extend to the matrix case. The basic result is: 

Theorem 5.5. The measures dvn o.nd Tr(Xn(x, x)) dfi{x) have the same weak limits. In 

particular, if djj, is regular, 

1 



lh:{Kn{x,x)dp,{x)) ^ dpE. (5.10) 
(n + 1)/ 

As in |174| . {-KuMxi^ny and (Tr^Mi-vr^) have a difference of traces which is bounded as n ^ oo, 
and this implies the result. Once one has this, combining it with Theorem 12.201 leads to: 

Theorem 5.6. Let I = [a,b] d E (ZM. with E compact. Let acss{dp) = E for anlxl matrix-valued 
measure, and suppose dfi is regular for E and 

dfi = W{x) dx + d/is (5.11) 

where dfig is singular and det{W{x)) > for a.e. x £ I. Then, 

(1) lim [p^{x)U^is{x)p^{x)^0, (5.12) 

n^oo J J 



(2) 



/ ;7TT T.Pfi^)^W{x)p,{x) - pe{x)1 



dx 0. (5.13) 
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5.5 Widom's Theorem 

Lemma 5.7. Let dfi be an I x I matrix-valued measure supported on a compact E <ZM. and let dr] 
be a scalar measure on R so 

dn{x) = W{x) dr]{x) + dfi^ix) (5.14) 
where dfig is dr] singular. Suppose for drj-a.e. x, 

dei{W{x)) > 0. (5.15) 

Then for drj-a.e. x, there is a positive real function C{x) so that 

||p^(x)|| <C(x)(n + 1)1. (5.16) 

In particular, 

\dei{p^{x))\<C{x)\n + l)'. (5.17) 

Proof. Since = 1, we have 

oo 

5^(n + l)-2||^,«|||<oo (5.18) 

n=0 

SO 

oo 

Y,{n + 1)-' Tt{p^{x)^W{x)p^Xx)) < oo 

for dr]-a.e. x. Since (j5.15p holds, for a.e. x, 

W{x) > b{x)l 
for some scalar function 6(x). Thus, for a.e. x, 



£(n+l)-2Tr(p^(x)tp^(x))<C( 



n=0 



Since < Tr(At^), we find (fSTTB]) . which in turn implies ([STTj) . □ 

This lemma replaces Lemma 4.1 of |171| and then the proof there of Theorem 1.12 extends to 
give (a matrix version of the theorem of Widom [195j): 

Theorem 5.8. Let dji be an I x I matrix-valued TTICQjSUTG with CXcss 

(dfi) = E cM. compact. Suppose 
dfi{x) = W{x) dpE{x) + djis{x) (5.19) 
with dfis singular with respect to dpE- Suppose for dpE-a.e. x, det(VF(x)) > 0. Then p is regular. 
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5.6 A Conjecture 

We end our discussion of regular MOPRL with a conjecture — an analog of a theorem of Stahl-Totik 
|180j : see also Theorem 1.13 of [171| for a proof and references. We expect the key will be some 
kind of matrix Remez inequality. For direct sums, this conjecture follows from Theorem 1.13 of 

m- 

Conjecture 5.9. Let -E be a finite union of disjoint closed intervals in M. Suppose |U is an / x / 
matrix-valued measure on M with. (Tggg (d/i) = E. For each t/ > and m = 1, 2, . . . , define 

Sm,, = {x : - + i]) > e-''"^!}. (5.20) 

Suppose that for each rj (with |-| = Lebesgue measure) 

lim \E \ SmJ = 0. 

m^oo 

Then fi is regular. 
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